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SECTION  I 
INTRODUCTION 

This  final  report  (in  two  volumes)  documents  the  results  of 
the  computational  fluid  dynamic  (CFD)  research  performed  by  the 
University  of  Dayton  for  the  Air  Force  Wright  Aeronautical 
Laboratories,  Aero  Propulsion  Laboratory  (AFWAL/PO) .  The  focus 
of  this  CFD  research  Is  the  confined,  turbulent  recirculating 
flowfield  behind  a  bluff  body  in  the  AFWAL/POSF  research 
combustor . 

In  an  existing  version  of  the  POSF  combustor,  a  centerbody 
configuration,  involving  the  turbulent  mixing  and  combustion  of 
an  annular  air  stream  and  a  central  fuel  jet  in  the  near  wake 
region  of  an  axisymmetric  bluff  body,  has  been  the  subject  of 
extensive  diagnostic  and  predictive  research.  The  results  of 
computations  of  this  configuration  for  the  nonreacting  flow  due 
to  the  annular  air  alone,  obtained  with  the  time-dependent 
Navier-Stokes  equations  incorporating  a  realistic  turbulence 
model,  are  presented  in  Volume  I  of  this  report.  The 
time-averaged  predictions  based  upon  a  solution  of  the 
Reynolds-averaged  Navier-Stokes  equations  for  a  proposed  POSF 
combustor  configuration,  Involving  two  annular  air  streams — a 
swirl-free  outer  stream  and  a  swirling  inner  stream,  a  central 
fuel  jet,  and  a  centerbody  imbedding  the  swirler  and  the  fuel 
nozzle,  are  reported  here,  in  Volume  II. 

1 .  BACKGROUND 

Research  in  gas  turbine  combustor-type  flows  is  motivated  by 
the  continuing  interest  in  enhancing  our  understanding  of  the 
turbulent  mixing  and  combustion  processes  in  turbojet-combustion 
chambers.  Thus,  the  POSF  research  combustor  (Reference  1)  has 
been  the  focus  of  extensive,  ongoing,  experimental  diagnostic  and 
computational  fluid  dynamic  research  activities.  These  activities 


have  been  directed  at  a  ducted  bluff-body  configuration  involving 
the  turbulent  flow  and  combustion  of  an  annular  air  stream  and  a 
central  fuel  jet.  Of  particular  interest  to  an  assessment  of 
time-dependent  calculations  for  combustor- type  flows  is  the  CFD 
research  completed  at  the  University  of  Dayton  Research  Institute 
(UDRI)  which  addressed  the  solution  of  time-averaged  and  time- 
dependent  Navier-Stokes  equations  for  the  nonreacting  flowfields 
in  the  POSF  centerbody  combustor  (References  2-3) .  A  further 
continuation  of  such  an  assessment  encompassing  time-dependent, 
turbulent  flow  computations  for  the  existing  centerbody  con¬ 
figuration  formed  part  of  the  present  Scholarly  Research  program. 
The  results  of  this  CFD  research  are  documented  in  Volume  I  of 
this  report. 

The  experimental  diagnostic  and  computational  predictive 
research  of  the  recent  past  on  the  existing  centerbody  configur¬ 
ation  has  given  rise  to  several  fundamental  and  practical 
implications  of  relevance  to  flowfields  in  gas-turbine  combus¬ 
tors.  A  practical  consequence  of  these  implications  is  an 
interesting  POSF  configuration  under  consideration  for  extensive 
development  and  testing.  The  new  POSF  configuration  is  expected 
to  simulate  more  closely  some  of  the  essential  flowfleld  features 
of  the  primary  zone  of  a  gas-turbine  combustor  than  did  the 
previous  centerbody  configuration. 

Although  the  proposed  POSF  combustor  continues  to  represent  a 
ducted  bluff-body  configuration,  it  is  distinguished  from  its 
predecessor  in  several  respects.  A  particularly  salient  feature 
is  the  presence  of  a  swirler  imbedded  in  the  bluff-body  flame 
holder  (which  also  carries  the  central  fuel  nozzle  as  before) . 
This  gives  rise  to  the  turbulent  mixing  of  two  annular  air 
streams — a  swirl-free  outer  stream  and  a  swirling  inner  stream. 
The  introduction  of  swirl  and  two  outer  streams  in  a  configu¬ 
ration  already  characterized  by  the  confining  outer  duct, 
bluff-body  recirculating  near  wake,  and  central  fuel  jet  renders 
the  new  combustor  extremely  complex.  Nevertheless,  it  imparts 
certain  realism  to  the  POSF  research  combustor. 


Accordingly,  a  CFD  examination  of  the  proposed  configuration 
was  considered  worthwhile  as  part  of  the  present  assessment 
studies.  Indeed,  such  an  investigation  conducted  prior  to  and 
during  the  POSF  development  program  could  provide  some  guidelines 
in  the  selection  of  optimum  conditions  for  further  development. 
Furthermore,  testing  of  CFD  methods  against  a  flowfield  that  did 
not  exist  and  had  not  been  measured  before  makes  it  possible  to 
ascertain  (when  new  measurements  become  available  for  comparison) 
whether  these  methods  could  really  serve  as  predictions  rather 
than  mere  "postdictions"  or  correlations. 

2.  SCOPE  OF  PRESENT  WORK 


As  part  of  an  assessment  of  time-dependent  calculations  for 
gas  turbine  combustor-type  flows,  a  proper  CFD  examination  of  the 
proposed  POSF  combustor  must  entail  both  unsteady  and  steady- 
state  flow  computations.  However,  within  the  limited  scope  of 
the  present  research  program  and  in  view  of  the  inadequate 
capability  of  the  time-dependent  formulation  (e.g.,  see  Volume  I) 
to  address  the  full  complexity  of  the  new  configuration,  only 
steady-state  calculations  were  carried  out.  Such  an  approach  is 
of  interest  nonetheless  for  this  configuration,  since  the 
numerical  solution  of  time-averaged  equations  appears  to  provide 
reasonable  predictions  of  the  overall  trends  in  gas  turbine 
combustor-type  flows  and  represents  an  accepted  design  tool  in 
the  industry  (Reference  4).  Moreover,  this  afforded  the  means 
for  a  computational  assessment  of  a  recent  time-averaged 
calculation  procedure  made  available  to  UDRI  by  the  Air  Force.  A 
newer  version  of  the  TEACH  procedure  which  seeks  to  solve  the 
Reynolds-averaged  Navler-Stokes  equations,  this  computer  program 
(Reference  5)  was  not  expected  to  involve  major  program 
development  effort.  Also,  this  version  appeared  to  possess  some 
physical  and  numerical  modeling  refinements  which  could  lead  to 
improved  flowfield  predictions  in  comparison  to  the  earlier 
results  (e.g.,  see  References  2  and  3). 


Thus,  the  present  CFD  research  involves  the  modification  and 
adaptation  of  the  time-averaged  formulation  of  Reference  5  for 
the  analysis  of  subsonic,  swirling,  and  reacting  turbulent 
flowfields  of  the  proposed  POSF  research  combustor.  The  plethora 
of  geometric  and  fluid  dynamic  parameters  characterizing  this 
complex  configuration  has  necessitated  a  somewhat  limited 
parametric  investigation.  As  will  be  seen  subsequently,  selected 
test  cases  were  considered  with  an  intent  to  emphasize  the  role 
of  the  inner  annular  stream  (with  or  without  swirl)  in  the 
turbulent  mixing  and  combustion  characteristics  of  the  POSF 
configuration . 

3.  OUTLINE  OF  REPORT 

The  governing  equations  and  the  boundary  conditions  are 
discussed  briefly  in  Section  II.  That  section  also  highlights 
the  computational  case  studies  and  the  parametric  values  thereof. 
The  results  of  the  CFD  investigation  are  presented  and  discussed 
in  Section  III,  where  some  of  the  deficiencies  of  the  calculation 
procedure  noted  in  the  present  investigation  are  also  Indicated. 
Section  IV  outlines  the  conclusions  from  the  present  study  and 
offers  our  recommendations  for  future  research. 


SECTION  II 

NUMERICAL  COMPUTATIONS 


This  section  deals  with  the  numerical  modeling  of  nonreacting 
and  reacting  turbulent  flowflelds  in  the  proposed  POSF  combustor 
through  the  framework  of  Reynolds-averaged  Navier-Stokes 
equations.  Following  is  a  brief  discussion  of  the  theoretical 
and  computational  aspects  of  the  solution  procedure.  Further 
details  are  available  in  Reference  5. 

1.  POSF  CONFIGURATION 

The  flowfield  of  Interest  to  the  present  CFD  research  is  the 
axlsymmetric  combustor,  a  schematic  of  which  is  seen  in  Figure  1. 
This  combustor  is  under  development  at  AFWAL/POSF  as  a  research 
tool  for  gas-turbine  combustor  modeling  and  diagnostic 
instrumentation  development.  In  contrast  with  the  centerbody 
combustor  configuration  of  earlier  studies  (References  1-3),  the 
configuration  of  Figure  1  Involves  an  additional  swirling  inner 
air  stream  which  is  bounded  between  the  nonswlrllng  outer  air 
stream  and  the  central  fuel  Jet.  With  the  introduction  of 
swirling  stream,  this  configuration  more  closely  resembles  some 
essential  aspects  of  gas-turbine  combustor  flowflelds.  The 
present  study  is  directed  at  an  evaluation  of  the  Influence  of 
the  swirling  stream  on  the  turbulent  mixing  and  combustion 
characteristics  of  the  bluff -body  near-wake  region. 

2 .  GOVERNING  EQUATIONS 

The  domain  of  the  present  CFD  Investigation  is  the 
open-ended,  cylindrical  chamber  downstream  of  the  base  of  the 
bluff  body.  Because  of  the  recirculating  nature  of  the  flow  in 
the  immediate  vicinity  of  the  bluff  body,  the  mathematical 
formulation  of  the  steady-state,  Reynolds-averaged  Navier-Stokes 
equations  entails  a  fully  elliptic  system  of  equations.  For  the 
axlsymmetric  configuration  of  Interest  here,  the  governing 
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Figure  1.  Scheaatic  of  the  Proposed  Combustor. 


equations  are  written  in  the  cylindrical  polar  coordinates  in  the 
following  general  form  for  all  the  dependent  variables: 

r  [fe  *  k  <"''♦>  -  «(rr*  -  k  <rr*  If']  ■  %  •“> 

Here,  9  denotes  the  time-mean  value  of  any  dependent  variable. 

U,  V,  and  H  are  the  respective  time-mean  velocity  components  in 

the  axial  (x) ,  radial  (r),  and  azimuthal  (6)  directions.  Because 

of  the  assumption  of  axlsymmetry  (in  the  mean),  Equation  (1) 

exhibits  no  explicit  dependence  on  6  for  all  the  flow  variables  4» 

(even  for  the  case  of  nonzero  swirl  velocity  W) .  p  is  the 

density,  r.  the  effective  exchange  coefficient  for  the  transport 
v 

of  the  variable  9,  and  S.  the  source  term  for  9.  The  latter 

9 

includes  true  source  terms  (such  as  those  arising  from  chemical 
reactions) ,  as  well  as  the  terms  not  covered  by  the  four 
convective  and  diffusive  terms  inside  the  square  brackets  of 
Equation  (1).  The  dependent  variable  9  denotes  U,  V,  W,  the 
turbulent  kinetic  energy  k,  its  dissipation  rate  e,  fuel  mass 
fraction  mf,  mixture  fraction  f,  and  enthalpy  h.  Table  l  shows 
all  the  relevant  information  for  these  dependent  variables. 

in  Table  1  denotes  the  appropriate  value  of  the  effective 
Prandtl /Schmidt  number  for  each  9. 

Turbulence  closure  is  obtained  from  the  two-equation  model  of 
Launder  and  Spalding  (Reference  6)  for  prescribing  the  eddy 
viscosity  from  the  two  scalars  k  and  c  (See  Table  1).  In  the 
computer  code  of  Reference  5  the  equations  for  mf  and  f  describe 
the  combustion  of  gaseous  propane  with  air.  The  mixture  fraction 
is  defined  by  the  sum  of  the  unburned  and  burned  propane  as 


where  the  replacement  of  the  mass  fraction  of  burned  propane  by 
one-third  of  the  mass  fraction  of  C02  follows  from  the 
stoichiometric  combustion  reaction  of  propane  and  oxygen  and  the 
assumption  of  equal  molecular  weights  for  propane  and  carbon 
dioxide.  The  enthalpy  Is  defined  by 

h  -  m^He  +  c  T,  (3 

r  K  p  ,  m 

where  H0  is  the  heat  of  reaction,  T  the  temperature,  and  a  th 

K  p ,  IB 

mean  specific  heat  of  mixture.  The  latter  is  obtained  from 

_  =  E  m.c  .,  (4 

P.ffl  j  j 

with  the  summation  being  over  all  species.  The  specific  heats 
cp  j  of  the  individual  species  are  obtained  from  polynomial  fits 
of  the  form 

2  3  4 

c  =(C  +CT+CT+CT+CT  )R/M,,  (5 

P'3  1J  2j  3  j  4j  5j  3 


where  R  is  the  universal  gas  constant,  Mj  is  the  molecular 
weights  of  species,  and  the  polynomial  constants  through 
Cg j  can  be  obtained  from  standard  thermochemical  tables  (see 
Reference  5  for  the  tabulated  values  for  the  species  CgHg ,  02 
N2,  H20  and  C02). 


The  source  terms  S.  for  all  $  are  seen  in  Table  l.  The 
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source  terms  Su  and  Sy  contain  respectively  the  axial  and  radial 
gradients  of  time-mean  pressure  p  which  is  one  remaining  unknown 
variable.  The  evaluation  of  p  requires  a  special  procedure  in 
the  solution  scheme  which  Is  discussed  in  Paragraph  11.4.  The 
source  terms  and  S£  correspond  to  the  standard  k-e  model  and 
do  not  include  the  additional  contributions  from  the  large 
streamline  curvature  and  the  preferential  influence  of  normal 
stresses  (which  were  Incorporated  in  References  2-3  previously) . 


The  source  term  Smf  for  fuel  mass  fraction  is  obtained  from 
the  Magnussen-Hjertager  combustion  model  (Reference  7).  This 
model  relates  the  rate  of  combustion  of  the  fuel  to  the  rate  of 
dissipation  of  turbulent  eddies.  The  source  terms  and  for 
mixture  fraction  and  enthalpy  vanish  identically  (because  of 
their  definitions). 

Finally  in  Table  1,  the  subscripts  1  and  t  in  p  and  a denote 
the  laminar  and  turbulent  counterparts  respectively.  Table  1 
also  lists  the  constants  appearing  in  the  turbulence  model,  as 
well  as  the  term  G  representing  the  rate  of  production  of 
turbulent  kinetic  energy  by  the  mean  motion.  Also,  note  tha*. 
once  mf  and  f  are  determined  by  the  solution  of  the  corresponding 
transport  equations  (for  4  ■  mf  and  f)  in  Equation  (l),  the  mass 
fractions  of  other  species  like  02,  N2,  C02  and  HgO  can  be  easily 
obtained  from  the  algebraic  relationships  available  in  Reference  5. 

3 .  BOUNDARY  CONDITIONS 

The  solution  of  Equation  (1)  requires  the  specification  of 
appropriate  boundary  conditions  for  each  of  the  dependent 
variables.  Since  the  governing  formulation  is  elliptic,  the 
boundary  conditions  are  prescribed  on  all  the  boundaries  of  the 
computational  domain.  With  symmetry,  only  the  top  half  of  the 
combustor  configuration  of  Figure  1  needs  to  be  represented  by 
the  computational  domain.  Thus,  the  top  boundary  is  the  duct 
wall  and  the  bottom  boundary  is  the  axis  of  symmetry.  The  left 
boundary  denotes  the  inflow  boundary  and  consists  of  the  outer 
annular  air  inlet,  bluff-body  face,  the  inner  annular  (swirling) 
air  inlet,  bluff-body  face,  and  the  central  fuel  inlet.  The 
right  boundary  is  the  outflow  boundary,  the  location  of  which  is 
unknown  a  priori.  Its  specification  is  arbitrary  and  it  is 
essential  to  ensure  that  the  sensitivity  of  the  computed 
solutions  to  the  specification  of  this  boundary  location  is  not 
significant. 


The  boundary  conditions  employed  in  the  present  CFD 
investigation  are  shown  in  Table  2.  Except  for  the  additional 
azimuthal  velocity  component  W,  the  present  boundary  conditions 
are  identical  to  those  used  in  our  previous  study  (Reference  2). 
Before  proceeding  further,  some  general  observations  concerning 
the  conditions  in  Table  2  can  be  made.  Along  the  bottom 
boundary,  3$/3r  *  0  for  all  $  and  v  =  0  due  to  the  symmetry 
requirement.  Along  the  top  boundary  and  bluff-body  face  (both 
representing  a  rigid  impermeable  wall),  the  normal  velocity 
components  and  the  normal  gradients  vanish,  the  tangential 
velocity  components  employ  the  law  of  the  wall,  and  the 
wall-function  formulations  also  apply  to  k  and  e.  Thus,  along 
duct  wall  V  and  diiij/dr  vanish  and  U,  W,  k,  and  £  utilize  wall 
functions.  Along  the  centerbody  face  U  and  Sn^/Sx  vanish  and  V, 
W,  K,  and  c  utilize  wall  functions.  At  the  outflow  boundary, 
consideration  of  overall  mass  conservation  governs  U,  while 
3 <|» / 3x  =  0  for  all  other  variables.  This  condition  for  the 
vanishing  axial  derivatives  at  the  exit  boundary  is  invariably 
used  in  time-averaged  computations  of  subsonic  internal  flows 
(e.g.,  see  References  2-4).  It  is  strictly  valid  only  for  fully 
developed  turbulent  flows  and  its  application  to  a  location  where 
the  flow  is  unlikely  to  be  fully  developed  has  always  been  a 
source  of  concern  to  the  practitioners  in  CFD.  The  unfavorable 
impact  of  this  application  on  the  interior  flowfield  development 
remains  to  be  addressed.  For  a  discussion  concerning  this  and 
other  exit-boundary  conditions  for  the  time-dependent 
formulation,  see  Reference  3. 

At  the  three  inlets,  all  the  variables  are  specified.  Thus, 
the  specification  of  Uin  follows  from  the  desired  mass  flow, 
provided  the  radial  velocity  V  is  zero.  Note  that  V  is  typically 
close  to  zero  for  well-designed  inlets.  The  swirl  velocity  is 
given  by 


Win  =  SR  Uin' 
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TABLE  2 

BOUNDARY  CONDITIONS 
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where  SR  is  the  swirl  ratio  which  is  presumed  known.  The  inlet 
value  of  k  is  prescribed  as  a  given  fraction  of  the  inlet  mean 
kinetic  energy.  Note  that  in  the  absence  of  data  for  turbulence 
intensity  in  all  three  orthogonal  coordinate  directions,  the 
assumption  of  isotropy  is  necessarily  invoked  in  specifying  k  as 
a  fraction  of  The  design  of  inlet,  therefore,  becomes  all 

the  more  critical  to  ensure  that  extreme  anisotropy  is  not 
present  there.  The  inlet  specification  of  e  follows  from 

1 • 5  0*75 

e  «  k  C  / ( 0 • 09A ) ,  (7 

in  n 

where  A  is  the  specified  boundary- layer  thickness  and  is  a 
constant  of  the  turbulence  model  (see  Table  1).  Note  that  this 
specification  differs  from  those  of  References  2-4  which  use  an 
inlet  turbulence-length  scale  related  to  the  radius  for  the 
central  inlet  and  to  the  difference  of  radii  for  the  annular 
inlet . 

Finally,  in  Table  2,  includes  m^ ,  f  and  h.  Since  no 
CC>2  is  present  in  the  inlet  stream,  according  to  Equation  (2),  f 
assumes  the  same  inlet  value  as  mf  (i.e.,  1  at  the  central  inlet 
and  0  at  both  the  annular  air  inlets) .  The  specification  of 
inlet  values  of  T  satisfies  the  corresponding  requirement  on  h 
[see  Equation  (3)].  The  vanishing  (axial  and  radial)  derivative 
conditions  at  other  boundaries  require  no  explanation,  except  for 
the  fact  that  the  use  of  the  Neumann  condition  on  h  at  the  duct 
wall  and  centerbody  face  implies  the  assumption  of  an  adiabatic 
wall.  For  cooled  walls,  the  specification  of  constant  wall 
temperature  should  be  used  as  the  boundary  condition  on  h  at  the 
solid  surface.  Thus,  Table  2  must  be  interpreted  appropriately 
for  the  variable  h. 


COMPUTATIONAL  DETAILS 


4  . 


The  application  of  the  TEACH  procedure  in  the  numerical 
treatment  of  the  Reynolds-averaged  Navier-Stokes  equations 
described  in  Paragraph  II. 2  is  well  known.  The  details 
concerning  the  underlying  theory  and  the  computational  procedure 
are  available  in  Reference  2  (for  the  standard  version)  and 
Reference  5  (for  the  refinements),  as  well  as  in  the  bibliography 
cited  therein.  For  the  present  discussion,  it  is  sufficient  to 
note  that  the  solution  procedure  (Reference  5)  involves  the 
Bounded  Skew-Upwind  Differencing  (BSUD)  method  of  Raithby 
(Reference  8)  in  place  of  the  hybrid  upwind  differencing  scheme 
used  previously.  The  calculation  of  the  pressure  field  which 
required  the  SIMPLE  algorithm  in  Reference  2  is  based  upon  the 
Pressure-Implicit  Split  Operation  (PISO)  in  Reference  5.  PISO 
involves  an  iterative  predictor-corrector  scheme  in  which  the 
velocity  is  first  calculated  for  a  guessed  pressure  distribution 
from  the  momentum  equations;  then  the  pressure  field  is  corrected 
so  as  to  have  the  velocity  field  satisfy  the  overall  continuity 
equation,  and  the  cycle  is  repeated.  Reference  5  must  be 
consulted  for  the  details  of  the  finite-difference  procedure. 

a .  Swirl  Characteristics 

Since  the  highlight  of  the  present  study  is  the  inclusion 
of  swirl  in  the  inflow  and  the  computation  of  the  azimuthal 
velocity  field  in  an  otherwise  axisymmetric  geometry,  we  will 
discuss  briefly  some  aspects  of  the  swirling  flows  relevant  to 
our  calculations. 

The  key  requirement  of  a  gas-turbine  combustor  is  that 
the  flame  remain  lighted  over  a  wide  range  of  operating  conditions 
This  is  usually  satisfied  by  the  introduction  of  swirl  which 
causes  recirculation  in  the  core  region  and  thereby  produces 
strong  shear,  high  turbulence,  and  rapid  mixing.  Moreover, 
swirling  flows  provide  increased  residence  time  for  the  fuel  in 
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a  combustor  of  given  length.  The  improvements  in  the 
flame-holding  characteristics  of  the  recirculation  zone  with 
swirl  are,  of  course,  accompanied  by  increased  pressure  drop. 
Also,  the  swirling  flows  are  known  to  lead  to  instability  in  the 
frequency  range  of  100-500  Hz  range. 

Beer  and  Chigier  (Reference  9)  define  the  swirl  number  as 
a  measure  to  characterize  the  amount  of  rotation  imparted  to  the 
axial  flow.  Thus,  the  swirl  number 

SN  =  2  Gm/(DsGt>'  (8) 

where  D  is  the  outer  diameter  of  the  swirling  jet  and  G  and 
s  m 

denote  respectively  the  integrated  axial  flux  of  the  angular 
momentum  and  the  axial  momentum  flux.  These  can  be  written  as 

D  /2 

G  =  /  2nrWpUdr  (9) 

m  0 

and 

D  /2 

G  *  /  2nrpU*dr  .  (10) 

0 

The  nature  of  the  flowfield  is  determined  by  the  magnitude  of  the 
swirl  number.  Thus,  for  SN  <  0*4,  there  is  no  recirculation  in 
the  potential  core  of  the  jet  and  the  swirl  is  said  to  be  weak. 
For  0*4  *  S„  <  0*6,  the  swirl  is  moderate  and  there  is 
considerable  divergence  of  the  streamlines  but  without 
recirculation.  S„  >>  0*6  denotes  strong  swirl,  a  condition 
characteristic  of  swirlers  of  practical  interest.  Note  that  in 
Equation  (8)  the  axial  momentum  of  both  the  central  and  annular 
streams  is  taken  into  account. 

In  the  combustor  configuration  of  Figure  1  much  of  the 
axial  momentum  is  contained  in  the  annular  streams,  especially 
the  outer  nonswirling  flow.  The  main  purpose  of  the  latter  is 
to  create  a  toroidal  recirculation  zone  behind  the  bluff  body. 


15 


>  V  V  V " 


r 

> 


The  strength  of  the  recirculation  zone  depends  not  only  on  the 
axial  momentum  of  the  outer  annular  stream,  but  also  on  the 
blockage  ratio  D^d^,  where  Dc  and  are  the  diameters  of  the 
centerbody  and  the  confining  duct  respectively.  The  flow 
reversal  of  the  central  jet  is  generally  affected  by  both  the 
recirculation  zone  of  the  outer  annular  stream  and  the  swirl  of 
the  inner  stream.  Unfortunately,  no  correlation  is  available  for 
the  relationship  between  the  flow-reversal  characteristics  of  the 
central  jet  and  the  social  momentum  of  the  annular  stream,  as  well 
as  the  blockage  ratio.  Therefore,  the  present  computational 
experiment  makes  use  of  a  parameter  called  the  Flux  Ratio  FR, 
defined  as  the  ratio  of  the  social  momentum  of  the  outer  annular 
stream  to  the  combined  momentum  of  both  the  central  jet  and  the 
inner  annular  stream. 

b.  Computational  Domain 

The  finite-difference  grid  employed  for  the  computational 
experiment  on  the  proposed  combustor  is  shown  in  Figure  2.  In 
the  construction  of  this  grid  with  37  axial  nodes  and  36  radial 
nodes,  an  exponential  stretching  was  used  in  the  axial  direction 
to  yield  a  denser  grid  near  the  face  of  the  bluff  body.  In  the 
radial  direction,  an  arbitrary  choice  of  6,  5,  and  11  mesh  points 
corresponds  to  the  central  jet  radius  and  the  inlets  of  the  inner 
and  outer  annular  streams.  We  believe  that  the  computational 
mesh  of  Figure  2  should  be  adequate  to  resolve  the  initial 
development  of  the  flowfield  with  reasonable  accuracy.  A 
comprehensive  experiment,  however,  is  required  to  investigate  the 
grid  sensitivity  of  the  computed  solutions. 

c .  Computational  Parameters  and  Test  Cases 


Some  of  the  constants  needed  by  the  governing  equations 
and  by  the  turbulence  and  combustion  models  therein  are  shown  in 
Table  1.  The  other  constants  used  in  the  present  computations 
are:  M,  *  1 • 88xlO“3Ns/m2 ,  Palr  *  1*1609  kg/m3,  Twall  *  293°K  and 
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Figure  2.  Finite-Difference  Computational  Grid. 


(for  the  combustion  of  propane)  HR  *  5*008  x  10'J/kg.  The  various 
parameters  used  in  the  different  test  cases  are  presented  in 
Table  3.  Unless  stated  otherwise  in  this  table,  the  Inflow  conditions 
used  are  the  following:  Pressure  p  =  98,000  Pa,  boundary- layer 
thickness  A  =  0*01  m,  k/Uin  =  0*005,  temperature  of  the  central 
jet  =  400*K,  and  the  temperature  of  the  outer  and  inner  annular 
streams  =  293 #K.  Note  that  in  Table  3  the  outer  and  inner  air 
streams  are  referred  to  as  the  annular  and  middle  flows. 

The  different  test  cases  we  examined  are  characterized  by 
different  values  of  the  axial  velocities  of  the  three  streams  and 
by  the  different  ratios  of  swirl-to-axial  velocities.  Of  course, 
only  the  middle  stream  has  nonzero  swirl  velocity.  The  central 
jet  consists  of  C02  or  propane,  depending  on  whether  nonreactive 
or  reactive  flow  is  considered.  Cases  1-8  of  Table  3  correspond 
to  the  former  and  cases  9-10  correspond  to  the  latter. 

The  parameters  describing  the  swirl  characteristics  are 
presented  in  Table  4.  Thus,  the  swirl  numbers  SN  considered  in 
the  present  study  varied  from  0*0  to  0*85.  Also  seen  in  Table  4 
are  the  values  of  the  flux  ratio  PR  [discussed  in  Paragraph 
II. 4. a]  used  for  the  different  cases.  FR  varies  from  3*76  to 
142*63.  Thus,  the  values  of  FR  for  cases  7-10  (corresponding  to 
an  outer  annular  mass  flow  of  2  kg/s)  are  an  order  of  magnitude 
larger  than  the  values  for  cases  1-6  (corresponding  to  a  mass 
flow  «  1  kg/s) .  This  distinction  gives  rise  to  the  annular- 
stream  dominant  flowfleld  and  the  Interactions  therein  with  swirl 
and  combustion  which  are  discussed  in  Section  III.  From  Tables  3 
and  4  it  is  easy  to  recognize  the  very  large  parametric  variation 
possible  in  the  computational  experiment.  The  limited  subset  of 
10  cases  we  have  investigated  merely  provides  the  flavor  in 
analyzing  the  effect  of  swirl  on  the  flowfleld  in  the  proposed 
combustor . 


TABLE  4 

SWIRL  NUMBERS  AND  FLUX  RATIOS  FOR  DIFFERENT  CASES  RUN 


CASE 


N 


0.7191 


0.4449 


0.8500 


0.7191 


0.7275 


p 

R 

34 

.59 

21 

.40 

10 

.22 

8 

.65 

3 

.76 

3 

.76 

142 

.63 

142 

.63 

142 

.63 

SECTION  III 

RESULTS  AND  DISCUSSION 


The  CFD  results  for  the  proposed  combustor  (see  Figure  1) 
obtained  from  calculations  outlined  in  the  previous  section  are 
discussed  in  the  following  paragraphs.  The  predictions  for  the 
nonreacting  flowfields  are  presented  first  and  those  for  the 
reacting  flowfields  next.  The  results  for  both  flowfields  are 
presented  and  discussed  in  terms  of  velocity-vector  plots  (and 
temperature-  and  density-contour  plots  for  reacting  flows),  cen¬ 
terline  (axial)  distributions  of  selected  flowfield  variables, 
and  radial  distributions  (at  different  axial  stations)  of  these 
variables. 

1.  NONREACTING  FLOWFIELDS 

Cases  1  through  8  (see  Table  3)  are  concerned  with  the 
nonreacting  flowfields.  The  complex  nature  of  the  flowfield 
interactions  in  the  near-wake  region  is  clearly  indicated  by  the 
computed  velocity  field.  The  following  discussion  of  the  overall 
flow  structure  discerned  from  the  velocity-vector  plots  sets  the 
stage  for  the  subsequent  examination  of  the  axial  and  radial 
profiles. 

a .  Velocity-Vector  Plots 

The  eight  cases  of  the  computed  velocity  field  considered 
here  correspond  to  only  a  limited  parametric  variation  of  several 
items  of  interest.  To  discern  quickly  the  flowfield  interactions 
when  several  variables  are  present,  the  computational  strategy 
adopted  in  these  eight  cases  is  to  keep  some  parameters  constant 
and  vary  others.  In  the  first  four  cases  the  outer  annular  mass 
flow  rate  (m  )  was  kept  constant  at  1  kg/s,  thereby  facilitating 
a  comparison  of  the  effects  of  variation  in  the  mass  flow  rates 
of  the  swirling  stream  and  the  central  stream.  In  the  next  four 
cases,  the  effects  of  variation  in  the  outer  annular  mass  flow 
and  in  swirl  were  examined  by  keeping  the  mass  flow  rates  of  the 
other  two  streams  constant . 
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Cases  1  through  4  consider  a  value  of  unity  for 
the  swirl  ratio  SR  (which  Is  defined  as  the  ratio  of  the 
tangential  to  the  axial  velocity  at  the  inlet)  of  the  Inner 
annular  stream.  The  velocity-vector  plots  for  these  four  cases 
are  seen  in  Figure  3.  Comparisons  of  cases  1  and  2  (Figures  3a 
and  3b),  and  of  cases  3  and  4  (Figures  3c  and  3d)  show  the 
effect  of  Increasing  the  central- jet  mass  flow  rate  (m )  at  two 
different  swirling-stream  mass  flow  rates  (m ) .  Comparisons  of 
cases  1  and  3,  and  of  cases  2  and  4  reveal  the  effect  of 
Increasing  m  at  two  different  values  of  m_. 

S  C 

Figure  3a  shows  that  the  wake  is  essentially 
characterized  by  four  vortices.  The  primary  vortex,  oriented  in 
the  clockwise  direction,  is  centered  at  the  normalized  axial  and 
radial  coordinates  of  0*45  and  0*325  respectively  (the 
normalization  being  with  respect  to  the  centerbody  diameter) . 

Two  smaller  vortices  of  opposite  sense  are  located  in  the  wake 
region  which  is  bounded  by  the  centerbody  face  and  the  primary 
vortex  in  the  axial  direction  and  by  the  centerbody  edge  and  the 
periphery  of  the  swirling  jet  in  the  radial  direction.  The 
clockwise  rotation  of  the  outer  vortex  and  the  counterclockwise 
rotation  of  the  inner  vortex  are  consistent  with  the  requirements 
of  the  overall  flowfleld  kinematics.  The  swirling  jet  is  drawn 
in  the  clockwise  direction  by  the  primary  vortex.  Finally,  a 
complete  flow  reversal  of  the  central  jet  occurs  with  an 
of f centerline  counterclockwise  vortex.  This  gives  rise  to  a 
forward  stagnation  point  on  the  centerline.  A  rear  stagnation 
point,  of  course,  signifies  the  end  of  the  recirculation  zone. 

The  normalized  centerline  locations  of  the  forward  and  rear 
stagnation  points  are  0*6  and  1*1  respectively. 

The  swirl  number  SR  for  this  case  is  0*7191.  In 
view  of  this  large  a  value  (S  >  0*6),  we  anticipate  that 
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Figure  3(b).  Velocity- Vector  Plots 


Figure  3(c).  Velocity- Vector  Plots 


Figure  3(d).’  Velocity- Vector  Plots. 


the  central- jet  flow  reversal  may  arise  from  the  swirling  stream 
alone.  This  is  so  because  with  sufficiently  strong  swirl,  a 
low-pressure  region  is  created,  thereby  resulting  in  a 
recirculation  of  the  core  region.  The  overall  effect  of  the 
swirling  stream  would  then  appear  to  be  the  creation  of  better 
mixing  characteristics  in  the  wake  region,  with  relatively  much 
smaller  mass  flow  in  the  outer  stream. 

A  comparison  of  Figures  3a  and  3b  shows  that  by 
Increasing  the  central- jet  mass  flow  rate  from  6  to  12  kg/hr,  the 
primary  recirculation  region  is  penetrated  by  the  central  jet. 
Although  the  centerline  flow  reversal  has  disappeared  in  Figure 
3b,  the  location  of  the  primary  vortex  center  remains  relatively 
unchanged  from  that  of  the  previous  case.  The  swirl  number  here 
is  0*44.  This  suggests  that  a  relatively  moderate  swirl  combined 
with  the  moderate  outer  annular  mass  flow  (1  kg/s)  is  not 
sufficiently  strong  to  cause  the  center line- flow  reversal. 

Farther  downstream,  the  central  jet  does  grow  in  size  since  it 
entrains  fluid  from  the  surrounding  region. 

Doubling  the  mass  flow  rate  of  the  swirling  stream 
(as  compared  to  Figure  3a  of  case  1)  results  in  a  flowfield  which 
is  dominated  by  the  swirling  stream,  as  seen  in  Figure  3c.  Table 
4  indicates  that  the  swirl  number  for  this  is  0*85.  The  primary 
vortex  is  seen  to  collapse  radially  inwards,  with  the  vortex 
center  lying  within  the  swirling  stream  at  a  normalized  distance 
of  approximately  0*355  from  the  centerbody  face.  The  normalized 
locations  of  the  forward  and  rear  stagnation  points  are  now  0*30 
and  1*105  respectively,  indicating  thereby  that  the  Increase  in 

SN  from  0*7191  to  0*85  halves  the  distance  of  the  centerline 
forward  stagnation  point. 

In  case  4  the  central- jet  mass  flow  Increases  to 
12  kg/hr  (as  compared  to  6  in  case  3)  and  decreases  to  0*7191. 
Figure  3d  shows  that  the  forward  stagnation  point  has  moved 
farther  downstream  to  0*68.  It  appears  that  the  swirl  is  very 


effective  in  causing  the  central-jet  flow  reversal  and  in 
enhancing  the  mixing  thereof.  Furthermore,  despite  the  higher 
value  of  central- jet  mass  flow,  flow  reversal  in  the  present  case 
has  been  achieved  by  increasing  the  swirling  mass  flow  from  0*028 
to  0*056  kg/s  and  the  swirl  number  from  0*44  to  0*7191  (as 
compared  to  case  2 ) . 

( 2 )  Variations  in  Swirl  and  Outer  Annular  Stream 

Cases  5  through  8  consider  constant  values  of  n»c 
(=  5*76  kg/hr)  and  mg  (=  o*028  kg/s).  A  comparison  of  cases  5 
and  6  reveals  the  effect  of  swirl  on  the  central-jet  flow 
reversal  at  a  low  value  of  outer  annular  mass  flow  (mg  =  o*328  kg/s) 
Cases  7  and  8  provide  this  comparison  at  the  much  higher  value  of 

ma  (=  2  kg/s).  Note  that  cases  5  and  7  are  swirl-free  (SR  =  0), 
while  cases  6  and  8  have  a  nonzero  swirl  (SR  =  i). 

The  presentation  of  the  velocity  fields  for  cases  5 
and  6  as  a  composite  plot  in  Figure  4a  facilitates  the  examina¬ 
tion  of  the  effect  of  swirl  on  the  central- jet  flow  reversal. 

The  bottom  half  of  the  plot  corresponds  to  case  5  and  the  top 
half  to  case  6.  The  swirl-free  flowfleld  in  case  5  shows  that 
there  is  no  centerline-flow  reversal  at  all.  The  flowfield  is 
seen  to  be  characterized  by  the  presence  of  two  weak  counter¬ 
clockwise  (it  should  be  kept  in  mind  that  this  is  so  only  below 
the  centerline)  vortices  in  the  wake.  The  outer  one,  of  course, 
is  the  primary  vortex  due  to  the  (outer)  annular  flow  past  the 
centerbody.  The  inner  vortex  is  located  very  close  to  the 
centerbody  face  between  the  middle  and  central  streams.  The 
incipient  tendency  for  entrainment  by  both  the  middle  and  central 
streams  is  also  noticeable  in  the  velocity  plot.  The  absence 
therein  of  two  additional  vortices  of  clockwise  direction  (again, 
only  below  the  centerline)  necessitated  by  kinematic  grounds, 
however,  may  be  due  to  an  inadequate  spatial  resolution  in  the 
near-wall  regions. 


Figure  4(a).  Velocity-Vector  Plots  (Top:  Case  6,  Bottom:  Case  5). 


The  introduction  of  swirl  in  the  middle  stream 
dramatically  changes  the  flowfield,  as  seen  in  the  top  half  of 
Figure  4a.  The  central  jet  proceeds  downstream  only  about  one 
radius  of  the  centerbody  before  encountering  a  complete  flow 
reversal.  Thus,  centerline  forward  and  rear  stagnation  points 
occur  at  the  normalized  distances  of  0*512  and  1*03  respectively. 
Furthermore,  a  comparison  of  the  top  and  bottom  halves  in  Figure 
4a  reveals  that  the  flowfield  structure  near  the  outer  vortex 
(i.e.,  radially  outside  the  swirling  stream)  remains  essentially 
unaffected  by  the  introduction  of  swirl,  whereas  rather  dramatic 
changes  are  noted  in  the  region  between  the  central  and  swirling 
streams.  Thus,  the  predominant  effect  of  swirl  is  felt  on  the 
central  jet  and  its  vicinity,  as  we  had  anticipated  (see  the 
discussion  in  Paragraph  II. 4. a). 

The  computed  velocity  fields  for  cases  7  and  8  are 
seen  in  the  composite  plot  of  Figure  4b.  The  bottom  half 
corresponds  to  the  swirl-free  case  7  and  the  top  half  to  the 
swirling  case  8.  Because  of  the  much  higher  mass  flow  in  the 
outer  annulus,  the  near  wake  is  dominated  by  the  large  primary 
vortex  and  characterized  by  the  centerline-flow  reversal, 
irrespective  of  the  presence  of  swirl  or  the  lack  thereof. 

Indeed,  the  outer  annular  flow  is  sufficiently  strong  to  cause 
the  reversal  of  both  the  central  and  swirling  streams.  Thus,  the 
comparison  of  cases  7  and  8  reveals  not  the  effect  of  swirl  in 
producing  the  centerline-flow  reversal  (as  was  noted  in  the 
comparison  of  cases  5  and  6),  but  the  subtler  effect  of  the 
degree  to  which  the  presence  of  swirl  modifies  the  outer- 
annular-stream  dominance.  This  observation  becomes  clear  from  an 
inspection  of  the  locations  of  the  vortex  centers  and  the 
centerline  stagnation  points. 

In  the  absence  of  swirl,  the  centerline  forward  and 
rear  stagnation  points  are  located  at  normalized  distances  of 
0*465  and  1*455  respectively.  With  the  introduction  of  swirl, 
the  forward  stagnation  point  has  moved  closer  to  the  bluff  body 
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and  is  at  the  normalized  distance  of  0*375.  Note  that  the 
location  of  the  rear  stagnation  point  has  remained  essentially 
unaffected.  That  the  predominant  effect  of  swirl  is  felt  only  In 
the  near  field  becomes  apparent  from  some  other  observations  as 
well.  Consider  one  of  the  secondary  vortices  located  closer  to 
the  bluff  body  and  between  the  outer  and  inner  annular  streams. 
This  is  a  clockwise  vortex  in  the  bottom  half  and  a  counter¬ 
clockwise  one  in  the  top  half.  The  introduction  of  swirl  pushes 
this  vortex  even  closer  to  the  bluff  body,  as  evidenced  by  the 
near  halving  of  both  the  axial  coordinate  of  the  vortex  center 
and  the  overall  size  of  this  vortex  (defined,  for  example,  by  the 
region  between  the  bluff-body  face  and  the  nearly  radially 
outward  velocity  vectors).  Note  that  the  axial  locations  of  the 
radially  outward  turning  of  the  middle  stream  and  a  corresponding 
stagnation  point  thereof  are  also  more  than  halved  by  the 
presence  of  swirl  in  the  middle  stream.  Similar  observations 
hold  for  the  other  secondary  vortex  occurring  between  the  central 
and  middle  streams  and  close  to  the  centerline  forward  stagnation 
point.  Finally,  the  introduction  of  swirl  even  affects  the 
primary  vortex  (albeit  to  a  marginal  extent).  The  normalized 
axial  and  radial  coordinates  of  the  vortex  center  are  0*575  and 
0*295  without  swirl  and  0*485  and  0*36  with  swirl  respectively. 

In  other  words,  with  the  introduction  of  swirl  the  primary  vortex 
is  also  pushed  upstream  towards  the  bluff  body  and  its  vortex 
center  moves  radially  outwards  and  axially  inwards.  Thus,  the 
overall  effect  of  swirl  is  to  reduce  the  spatial  extent  of  the 
mixing  region  and  thereby  enhance  the  rate  of  mixing  therein. 

b.  Centerline  Variations 

Having  seen  the  overall  structure  of  the  P0SF  combustor 
flowfield  in  the  eight  cases,  we  turn  our  attention  to  the  axial 
variation  of  selected  flow  variables  along  the  centerline. 

Present  discussion  is  brief  and  the  prospective  reader  is 
encouraged  to  interpret  the  noted  flowfield  behavior  in  this 
discussion  from  an  integrated  perspective  derived  from  the 


overall  flowfield  description  presented  earlier  and  the  radial 
distributions  discussed  subsequently. 

( 1 )  Axial  Velocity  Fields 

The  predicted  centerline  variation  of  the  mean  and 
root  mean  square  (rms)  axial  velocity  components  for  the  eight 
cases  is  seen  in  Figures  5a  through  5h.  Both  the  mean  (U)  and 
the  rms  (u)  components  are  normalized  by  the  mean  inlet  axial 
velocity  (Ua)  of  the  outer  annular  stream.  Note  that  the 
assumption  of  isotropy  is  Invoked  in  deriving  u  as  given  by 
(2/Sk)1/2.  The  abscissa  in  the  plots  denotes  the  axial  distance 
(x)  normalized  by  the  centerbody  diameter  (D) . 

In  all  eight  cases,  the  mean  axial  velocity  remains 
constant  for  approximately  two  diameters  of  the  central  jet 
before  it  shows  a  linear  decay.  This  linear  dependence  is 
observed  for  2  «  x/Dc  <  6  to  7*5,  where  Dc  is  the  initial 
diameter  of  the  central  jet. 

For  a  jet  issuing  from  a  nozzle  into  the  ambient 
atmosphere,  the  initial  distance  over  which  the  mean  axial 
velocity  component  and  other  scalar  variables  remain  invariant  is 
the  potential  core.  For  the  present  configuration,  the  predicted 
potential  core  of  the  central  jet  extends  to  about  2  jet  diame¬ 
ters  when  both  the  outer  and  inner  annular  streams  are  present 
and  there  is  swirl  in  the  middle  stream  (this  is  also  true  of  the 
zero-swirl  cases  5  and  7).  In  the  previous  ducted  centerbody 
configuration  laser  Doppler  anemometry  (LDA)  measurements 
(Reference  10)  showed  that  the  potential  core  of  the  central  jet 
(in  the  presence  of  very  low  annular  air  flow  which  existed 
merely  to  carry  the  seed  particles)  extended  to  6 *9-7 *5  jet 
diameters  in  the  mass  flow  range  of  6  to  12  kg/hr.  Clearly,  in 
the  new  POSF  configuration  the  addition  of  two  outer  streams  and 
swirl  has  significantly  altered  the  entrainment  characteristics 
and  reduced  the  length  of  the  potential  core. 
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Figure  5(b),  Centerline  Profiles  of  Mean  and  raa  Axial  Velocity 
Coaponents . 
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To  the  extent  that  earlier  measurements  (Reference 
10)  confirmed  that  the  development  of  the  central  jet  issuing 
from  a  vertical  wall  (i.e.,  the  bluff-body  face)  in  a  confined 
situation  initially  resembled  free-jet  behavior,  the  present 
predictions  indicate  the  strong  influence  of  the  annular  streams 
on  the  initial  development  of  the  central  jet.  Let  us  denote  the 
present  potential-core  length  by  L  and  the  previous  free-jetlike 
potential-core  length  by  Lf.  What  is  interesting  to  note  in  the 
present  predictions  is  that  L^  is  very  close  to  the  sum  of  L  and 
the  present  distance  of  linear  decay.  The  reason  for  this 
coincidence  is  not  clear. 

Figures  5a  through  5h  further  show  that  the  center- 
line  location  signifying  the  end  of  the  linear  decay  of  the  mean 
velocity  appears  to  coincide  with  the  location  of  the  peak  rms 
velocity  (or  of  the  first  peak  in  cases  with  more  than  one  peak) . 
Farther  downstream  of  this  peak,  the  rms  velocity  starts  to  fall 
until  the  centerline  location  of  the  peak  negative  mean  velocity 
is  reached.  Beyond  this  minimum,  the  rms  velocity  Increases 
slightly  before  it  tapers  off  gradually  towards  the  exit 
boundary.  For  the  two  cases  (viz. ,  2  and  5)  where  the 
central-jet  penetration  has  eliminated  the  centerline-flow 
reversal,  the  rms  velocity  falls  gradually  towards  the  exit 
boundary  after  its  peak  at  Lf  (see  Figures  5b  and  5e) . 

The  two  cases  where  the  centerline  profile  of  rms 
velocity  exhibits  two  sharp  peaks  correspond  to  cases  7  and  8  of 
higher  mass  flow  (2  kg/s)  in  the  outer  annular  stream  (see 
Figures  5g  and  5h) .  While  the  first  peak  occurs  at  Lf,  the 
second  peak  occurs  very  slightly  downstream  of  the  forward 
stagnation  point.  The  latter  occurrence  is  certainly  reminiscent 
of  the  behavior  noted  in  both  predictions  (References  2-3)  and 
measurements  (Reference  10)  of  the  previous  POSF  centerbody 
configuration.  The  second  peak  in  the  present  predictions  is 
clearly  due  to  the  higher  momentum  flux  of  the  outer  annular 
stream.  However,  only  experimental  data  in  the  new  configuration 
can  quantitatively  validate  the  present  predicted  behavior. 
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For  cases  with  centerline-flow  reversal  (see 


Figures  5a,  5c,  5d,  5f,  5g,  and  5h) ,  the  mean  axial  velocity 
increases  downstream  of  the  location  of  the  peak  negative 
velocity  and  past  the  rear  stagnation  point.  The  quantitative 
accuracy  of  the  predicted  rate  of  recovery  beyond  the  rear 
stagnation  point,  however,  remains  to  be  verified  by 
measurements . 

( 2 )  Temperature  and  Concentration  Fields 

Figures  6a  through  6h  show  the  centerline  variations 

of  CC>2  mass  fraction  and  temperature.  Here,  the  temperature  is 

the  normalized  variable  given  by  (T  -  T  .  /(T  -  T  .  ), 

M  1  '  min  '  v  max  mm7 

T  being  the  maximum  value  of  temperature  in  the  flowfield. 
max 

For  the  present  results,  T  is  400“K,  corresponding  to  the 

nidX 

inflow  temperature  of  CC>2 . 

The  one  thing  readily  apparent  in  Figures  6a  through 
6h  is  that  with  the  above  normalization  of  the  temperature,  the 
centerline  profiles  of  both  temperature  and  CO 2  mass  fraction  are 
nearly  identical.  Of  course,  this  behavior  is  a  direct  conse¬ 
quence  of  the  assumption  of  unity  Lewis  numbers  for  both  laminar 
and  turbulent  transport  coefficients  of  the  scalar  fields  in  the 
numerical  computation  of  these  nonreacting  cases. 

An  interesting  aspect  of  the  present  predictions 
that  is  less  readily  apparent  relates  to  the  similarity  of  the 
mean  axial  velocity  profiles  of  Figure  5  and  the  corresponding 
scalar  profiles  of  Figure  6.  Had  the  mean  velocity  in  Figure  5 
been  normalized  with  respect  to  the  central-jet  inlet  velocity 
U  [instead  of  the  annular  velocity  U  as  was  done  in  Paragraph 
Ill.l.b(l)],  we  would  readily  see  that  the  mean  axial  velocity 
profile  also  is  nearly  identical  with  the  scalar  profiles  for 

x/p  ^  Lf.  Thus,  the  profiles  of  C09  mass  fraction,  normalized 
c  r  i 

temperature,  and  U/U  would  remain  constant  over  a  distance  of 

about  2D  from  the  centerbody  face  and  then  exhibit  a  linear 
c 

decay  for  about  6  to  7«5  D  .  This  linear  dependence  can  be 
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Figure  6(a) .  Centerline  Profiles  of  CO  Mass  Fraction  and 
Temperature .  2 
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Centerline  Profiles  of  CO  Mess  Fraction  and 
Teaperature .  2 


47 


>  V 
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quantified  by  the  following  correlation: 


m__  =  (T  -  T  .  )/ (T  -  T  ,  )  =  U/U  = 
C02  '  mm  max  min  c 


1  -  0*44(x/Dc-2)/(Lf-2) 


(11) 


It  would  be  interesting  to  see  if  the  proposed  measurements  in 
the  new  POSF  configuration  confirm  the  correctness  of  Equation 
(11)  • 

( 3 )  Implications  of  Present  Predictions 

Further  examination  of  the  implications  of  the 
present  predictions  of  the  axial-velocity  and  scalar  fields  is 
instructive.  That  the  development  of  the  central  jet  in  the 
proposed  configuration  departs  significantly  from  the  free-jet¬ 
like  behavior  noted  in  Reference  10  should  not  be  surprising.  To 
be  sure,  the  present  departure  stems  from  the  presence  of  the 
outer  and  inner  annular  streams,  as  well  as  of  swirl  in  the 
latter  (except  for  cases  5  and  7).  These  were  absent  under  the 
conditions  of  the  previous  measurements  (Reference  10)  which 
addressed  the  development  of  a  central  jet  in  a  negligibly  small 
coflowing  annular  stream  (present  solely  for  carrying  the  LDA 
seed  particles) .  Indeed,  earlier  CFD  predictions  of  the  central- 
jet  dominant  flowfields  also  indicated  free-jetlike  behavior 
(see,  e.g.,  References  2,  11,  and  12).  Note  that  the  central- 
jet  dominant  flowfields  in  the  proposed  POSF  configuration  also 
exhibit  free- jet  like  behavior.  This  becomes  clear  from  an 
inspection  of  Figures  5  and  6  (b  and  e)  corresponding  to  cases  2 
and  5  where  the  central  jet  penetrates  the  recirculation  zone  and 
eliminates  the  centerline-flow  reversal.  The  profiles  for  the 
mean  axial  velocity,  C02  mass  fraction,  and  temperature  exhibit 
the  characteristic  1/x  -  dependence  for  a  distance  of  about  30 
initial  diameters  of  the  central  jet.  Presumably,  the  free-jet 

behavior  becomes  the  better,  the  smaller  the  values  of  m  and  m  . 

a  s 
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Until  experimental  data  from  the  proposed  combustor 
are  available,  we  cannot  ascertain  the  validity  of  the  myriad 
implications  from  the  present  predictions.  Furthermore,  the 
computational  procedure  of  Reference  5  itself  has  not  been  tested 
as  comprehensively  as  the  earlier  procedures  were  in  References 
2-4  and  12.  Recently  completed  CFD  research  (Reference  13)  did 
attempt  a  comparative  study  of  the  predictions  from  the  old  and 
new  procedures,  with  respect  to  the  measurements  in  the  old 
centerbody  combustor .  However ,  that  study  only  addressed  a 
limited  inquiry  concerning  the  primary  vortex  center  and 
concluded  that  the  new  computational  procedure  was  superior  in 
its  predictions.  Therefore,  it  is  of  interest  to  further 
scrutinize  the  present  predictions.  It  must  be  stressed, 
however,  that  the  discussion  here  must  be  regarded  as  speculative 
at  best.  Note  that  for  the  proposed  configuration  experimental 
data  are  not  available,  and  predictions  have  not  been  made  with 
the  earlier  computational  procedure.  So  what  is  attempted  here 
is  a  comparative  study,  the  conclusions  of  which  may  turn  out  to 
have  dubious  validity. 

(a)  Aspects  of  Conformity 

For  this  scrutiny,  we  consider  the  swirl-free 
case  7  which  has  an  outer  annular  mass  flow  of  2*01  kg/s,  and  a 
central- jet  mass  flow  of  5*76  kg/hr.  Although  swirl-free,  the 
inner  annular  stream  has  the  nonzero  mass  flow  rate  of  0*028  kg/s. 
We  compare  the  present  prediction  for  this  case  with  the 
prediction  of  the  standard  k-e  model  for  2  kg/s  annular  air  flow 
and  6  kg/hr  C02  central  flow  in  Reference  2.  Figure  21  of 
Reference  2  compares  the  prediction  and  measurement  (of  Reference 
10)  for  the  centerline  profiles  of  the  mean  and  rms  axial 
velocity  components.  The  normalizations  of  the  abscissa  and  the 
ordinate  in  this  figure  are  the  same  as  in  Figure  5g.  It  is 
noteworthy  that  the  overall  trends  of  the  mean  and  rms  velocity 
profiles  in  the  old  and  new  POSF  configurations  are  similar.  The 
one  distinction  pertains  to  the  rms  velocity  profile:  the 


prediction  for  the  new  configuration  has  two  sharp  peaks,  whereas 
there  was  only  one  in  the  earlier  prediction  (l.e.,  old  configu¬ 
ration)  . 

There  is  reason  to  believe  that  the  mean 
velocity  profile  in  Figure  5g  will  show  good  agreement  with  the 
proposed  measurement.  Although  the  initial  velocity  ratio  at 
present  is  1*4  (in  contrast  with  the  value  of  1*15  in  Reference 
2),  this  is  a  consequence  of  the  different  Inflow  areas  in  the 
new  configuration.  This  change  in  the  velocity  ratio,  of  course, 
explains  the  Increase  in  the  normalized  location  of  the  forward 
stagnation  point  from  about  0*28  (measurement)  earlier  to  about 
0*46  (prediction)  now  (note  that  there  is  no  linear  scaling, 
however).  The  Increase  in  the  normalized  location  of  the  rear 
stagnation  point  from  about  1  (both  measurement  and  prediction) 
earlier  to  about  1*44  (prediction)  now  is  not  easily  explained  by 
the  increase  in  the  velocity  ratio.  The  presence  of  the  middle 
stream  could  conceivably  have  a  significant  influence  in  this 
regard.  Note,  however,  the  ratio  of  the  rear-to-forward 
stagnation  point  distances  in  the  proposed  configuration  is  3*13 
(prediction) .  That  this  is  fairly  consistent  with  the  value  of 
3*57  (measurement)  in  the  old  configuration  provides  some 
quantitative  basis  for  our  hope  that  the  present  mean-velocity 
prediction  will  show  reasonable  agreement  with  the  proposed 
measurement . 

A  comparison  of  the  centerline  CC>2  mass 
fraction  profiles  in  the  old  and  new  configurations  (Figure  23 
from  Reference  2  and  Figure  6g)  again  brings  out  the  similarity 
in  the  overall  trends.  As  in  the  case  of  the  axial  mean  velocity 
profiles  (and  certainly  because  of  the  differences  therein),  there 
are  quantitative  differences  between  the  old  [measurement 
(Reference  14)  and  prediction]  and  the  new  (prediction) 
configurations.  And  some  of  these  differences  could  be  real. 
Nevertheless,  further  insights  can  be  gleaned  from  the  scalar 
profiles.  Reference  2  noted  how  its  predictions  and  the 


experimental  measurements  (References  10  and  14)  confirmed  the 
earlier  anticipated  trends  (Reference  11),  of  a  rapid  decay  first 
and  an  equally  rapid  approach  subsequently  to  uniform  values, 
exhibited  by  the  centerline  CO^  concentration.  Reference  2  had 
also  noted  that  the  point  of  intersection  obtained  by  extrapo¬ 
lating  the  portions  of  the  profiles  that  denote  the  rapid  decay 
and  the  approach  to  unif ormizat ion  in  Figure  23  fell  very  close 
to  the  forward  stagnation  point  in  Figure  21.  Indeed,  a 


rationale  for  expecting  precisely  this  occurrence  had  been 
suggested  in  Reference  11.  It  is  appealing  that  the  present 
predictions  also  display  this  behavior  for  the  annular-stream- 
dominant  regime.  An  inspection  of  Figures  5g  and  6g  (the  linear 
segment  is  extrapolated  beyond  to  intersect  with  the  extrapo¬ 
lation  of  the  nearly  horizontal  segment)  reveals  this  internal 
consistency  of  the  present  predictions.  Thus,  we  conclude  that 
the  present  predicted  behavior  of  the  mean  axial  velocity,  CC>2 
mass  fraction,  and  temperature  fields  should  show  reasonable 
agreement  with  the  yet-to-be-available  measured  results.  An 
experimental  verification  of  the  predicted  behavior  will  be 
gratifying  in  view  of  the  "true"  predictions  of  the  present 
research . 


(b)  Area  of  Disagreement 


It  is  very  unlikely  that  the  foregoing  will 
extend  to  the  rms  velocity  field.  Figure  21  of  Reference  2  has 
clearly  shown  that  the  predicted  rms  axial  velocity  profile 
compares  rather  poorly  with  the  LDA  results  of  Reference  10. 

This  discrepancy  between  the  prediction  and  measurement  is 
especially  significant  immediately  downstream  of  the  bluff  body. 
The  assumption  of  isotropy  in  the  calculations  in  the  face  of 
considerable  anisotropic  effects  in  the  near  wake  is  known  to 
contribute,  at  least  in  part,  to  the  observed  discrepancy.  There 
are  other  aspects  of  weakness  in  the  standard  k-e  model,  such  as 
the  nonaccounting  for  the  extra  strain  rates  inherent  in  the 
large  streamline  curvature  of  the  flowfield,  the  noninclusion  of 
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the  preferential  influence  of  the  normal  stresses  in  the  e- 
equation,  etc.  It  should  be  noted,  however,  that  as  seen  in 
Figure  22  of  Reference  3,  our  ad  hoc  corrections  for  these 
effects  did  not  improve  the  nature  of  the  rms  velocity 
predictions . 

Thus,  we  believe  that  the  rms  velocity  profile 
predicted  in  Figure  5g  may  not  show  as  good  an  agreement  with  the 
experimental  profile  as  we  anticipate  for  the  mean  field  compar¬ 
isons.  Such  an  outcome  need  not  cause  surprise  or  concern.  All 
the  earlier  CFD  research  by  us  and  others  to  address  complex 
recirculating  turbulent  flowfields  by  means  of  the  Reynolds- 
averaged  formulation  (employing  any  eddy  viscosity  derived  from 
the  k-e  model)  has  been  confronted  by  this  dichotomy  of 
reasonable  prediction  of  the  mean  field  and  poor  prediction  of 
the  fluctuating  field.  Therefore,  it  is  worthwhile  to  ask  why 
the  mean-field  predictions  are  reasonable  and,  in  fact,  are 
better  than  could  be  justifiably  expected  from  the  k-e  model. 

The  answer  to  this  question  could  be  given  at  two  different 
levels . 

(c)  Rationale 

On  a  practical  computational  level,  we  could 
accept  the  position  (advanced  by  Peter  Bradshaw  of  the  Imperial 
College,  among  others)  that  the  turbulent-flow  behavior  is  so 
complex  and  the  search  for  general  closures  is  so  futile  that 
the  CFD  study  of  realistic  flows  will  require  the  development  of 
geometry-specific  empiricism.  Indeed,  this  viewpoint  has  been 
voiced  earlier  by  Reference  15  which  invoked  Bradshaw's  arguments 
that  in  complex  flows  driven  essentially  by  static  pressure 
gradients,  the  turbulent  shear  stresses  are  at  least  two  orders 
of  magnitude  smaller  than  typical  dynamic  pressures.  Thus,  while 
the  details  of  turbulent  stress  gradients  are  locally  important 
and  their  accurate  predictions  will  require  correct  modeling 
assumptions,  the  global  flowfield  is  relatively  insensitive  to 


model  inaccuracies .  Furthermore,  as  emphasized  in  Reference  15, 
confined  flowfields  are  characterized  by  large  turbulence 
production  rates  at  the  solid  boundaries  and  by  the  presence  of 
small  length  scales  (near  these  boundaries)  which  provide  a  ready 
sink  for  the  turbulence  energy.  This  fact  generally  leads  to  a 
near  equilibrium  between  the  production  and  dissipation  of 
turbulence  energy  over  most  of  the  flowfield  and  ensures  the 
"success"  of  the  k-e  model  in  confined  flows  (Reference  16). 


On  a  more  fundamental  level,  the  probable 
explanation  of  how  much  can  be  computed,  and  how  well,  using 
relatively  simple  assumptions  regarding  the  turbulence  is  offered 
by  Lumley  (Reference  17).  According  to  Lumley,  the  plausibility 
of  a  turbulence  model,  even  while  representing  not  quite  real 
turbulence,  is  ensured  if  it  conserves  momentum  and  energy; 
transports  the  right  amount  of  everything  (e.g.,  momentum, 
energy,  Reynolds  stress,  heat  flux,  etc.),  although  not  by  quite 
the  right  mechanism;  satisfies  thermodynamic  realizability  [in 
the  sense  of  Schumann  (Reference  18),  according  to  which 
nonnegative  quantities  are  never  negative,  Schwarz'  inequality  is 
always  satisfied,  etc.];  behaves  correctly  for  both  large  and 
small  Reynolds  numbers;  and  reduces  to  real  turbulence  in  one 
limit  (weak  inhomogeneity  and  unsteadiness) .  Presumably,  any 
model  that  satisfied  all  these  restrictions  would  behave  about 
the  same  and  all  that  remains  is  to  fix  the  amount  of  transport 
through  physical  input  from  experiment. 

Thus,  it  would  seem  that  in  confined 
flowfields,  such  as  those  exemplified  by  the  old  and  new  POSF 
combustor  configurations,  the  turbulence  is  simply  a  mechanism 
for  momentum  transfer,  and  usually  downhill.  Modeling  the 
turbulence  by  any  reasonable  mechanism  (such  as  the  Reynolds- 
averaged  formulation  and  an  eddy-viscosity  model)  which 
transports  about  the  right  amount,  and  which  guarantees  thermo¬ 
dynamic  realizability  will  be  adequate  to  describe  the  global 
flowfield,  so  long  as  the  mean  motion  is  handled  correctly.  To 
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be  sure,  this  simplistic  approach  has  not  been,  and  will  not  be, 
successful  in  predicting  the  finer  details  and  where  additional 
complexities  such  as  chemical  reactions  occur. 


(4)  Swirl  Velocity  Fields 

The  axial  variation  of  the  mean  swirl  velocity 
component  for  the  six  nonzero-swirl  cases  (i.e.,  all  cases  except 
cases  5  and  7)  considered  is  seen  in  Figures  7a  through  7f. 

Since  the  mean  swirl  velocity  component  should  vanish  on  the 
centerline,  the  axial  variation  seen  in  the  figures  is  not  the 
centerline  profile  but  represents  the  computed  swirl  velocity 
close  to  the  centerline.  The  swirl  velocity  is  shown  normalized 
with  respect  to  the  inflow  swirl  velocity  (W  )  of  the  inner 
annular  stream.  The  plots  also  show  the  previously  discussed 
centerline  variation  of  COg  ma8S  fraction  [Paragraph  III.l.b(2)] 
as  a  frame  of  reference  for  the  discussion  of  the  swirl-velocity 
profiles. 


The  rather  small  magnitudes  of  the  normalized  swirl 
velocity  components  (<  1*  except  for  case  2)  in  the  vicinity  of 
the  centerline  should  be  clear  from  the  figures.  The  use  of  a 
scale  factor  of  10  or  100  in  the  ordinates  helps  one  to  discern 
the  significant  axial  variations  in  the  normalized  swirl 
velocity.  From  the  viewpoint  of  experimental  verification,  the 
vicinity  of  the  centerline  may  not  be  the  optimum  location  for 
examining  the  axial  variation  of  the  mean  swirl  velocity.  As 
will  be  seen  subsequently  in  our  discussion  of  the  radial 
distributions  of  the  swirl  velocity  [see  Paragraph  III.l.c(3)], 
the  predicted  profile  in  the  axial  direction  2*5  cm  radially 
outward  from  the  centerline  will  be  more  helpful  in  experimental 
verification.  At  this  location  (which  corresponds  to  the  radial 
coordinate  of  the  swirler  axis)  the  magnitudes  of  the  swirl 
velocity  components  are  much  larger.  The  present  discussion, 
however,  does  not  include  such  an  of f centerline  profile. 
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It  is  clear  that  for  all  six  cases  the  normalized 
swirl  velocity  first  peaks  just  downstream  of  Lf .  Subsequently, 
it  falls  off  sharply  to  a  minimum  at  a  location  corresponding 
to  the  forward  stagnation  point  (except  in  Figure  7b  of  case  2 
where  the  central  jet  penetrates  the  recirculation  region  and 
eliminates  the  centerline-flow  reversal).  There  is  a  sharp  rise 
to  a  second  peak  in  the  swirl  velocity  beyond  the  minimum. 

Farther  downstream  there  is  a  gradual  tapering  off  towards  the 
exit  boundary  in  all  the  cases.  The  magnitude  of  the  second  peak 
relative  to  the  first  peak,  however,  appears  to  depend  on  other 
parameters  of  the  problem  (such  as  m  ,  m  .  m _,  SM,  etc.).  Thus, 
in  cases  1,  4,  6,  and  8  the  first  peak  represents  the  absolute 
maximum,  while  case  3  (see  Figure  7c),  corresponding  to  the 
maximum  swirl  number  0*85  considered  in  this  study,  exhibits  the 
second  peak  as  the  absolute  maximum  (nearly  three  times  as  large 
as  the  first  peak) . 

A  comparison  of  Figures  7a  and  7c  shows  that  the  near 

wake  is  dominated  by  the  swirling  stream  when  its  mass  flow  rate 

m  is  doubled.  This  is  evident  both  from  the  shift  towards  the 
s 

bluff-body  face  in  the  locations  of  the  pea^s  and  from  the 

decrease  (by  more  than  a  factor  of  2)  in  the  magnitude  of  the 

first  peak  and  the  increase  (by  more  than  a  factor  of  4)  in  the 

magnitude  of  the  second  peak.  An  increase  in  the  central  jet 

mass  flow,  with  the  mass  flow  rates  in  the  two  outer  streams  kept 

fixed,  diminishes  the  magnitudes  of  both  peaks  in  swirl  velocity 

considerably.  This  is  clear  from  a  comparison  of  Figures  7c  and 

7d  which  indicates  that  the  increase  in  m  is  to  counter  the 

c 

influence  of  the  swirling  stream.  Qualitatively,  the  trends  in 
Figures  7a  and  7d  have  become  similar.  This  is  consistent  with 
the  identical  value  of  swirl  number  (see  Table  4)  in  cases  1  and 
4  . 

One  may  expect  that  a  decrease  in  the  outer  annular 
mass  flow  rate  (m  )  would  produce  a  result  somewhat  similar  to 

a 

that  of  an  increase  in  m  .  This  is  the  case,  as  shown  by  a 

c 


shown  by  a  comparison  of  Figures  7a  and  7e.  There  is  a 
significant  reduction  in  the  magnitudes  of  both  peaks.  The 
relative  increase  in  the  second  peak  vis-a-vis  the  first  peak  is 
larger,  however.  Finally,  an  increase  in  m  from  1  to  2  kg/s 
essentially  preserves  the  trend.  This  is  seen  from  a  comparison 
of  Figures  7a  and  7f  which  indicates  that  the  first  peak  remains 
unchanged  (both  in  magnitude  and  location) ,  whereas  the  second 
peak  is  decreased  in  magnitude.  It  is  worth  noting  that  the 
similarity  of  the  trends  in  cases  1,  4,  6,  and  8  is  consistent 
with  the  nearly  identical  values  of  SN< 

c.  Radial  Distributions 


A  further  examination  of  the  predicted  behavior  of  the 
POSF  combustor  flowfields  can  be  made  by  considering  the  radial 
variation  of  selected  flow  variables.  For  this  purpose  four 
axial  stations  were  selected  (located  at  0*1,  0»73,  1*58,  and 
4*12  cm  from  the  centerbody  face)  and  the  radial  distributions 
thereat  of  the  turbulent  kinetic  energy,  temperature,  C02  mass 
fraction,  and  mean  swirl  velocity  were  studied.  In  the  plots  of 
the  profiles  both  the  axial  and  radial  distances  are  normalized 
with  respect  to  the  centerbody  radius  H  (=  7  cm) .  The  normalized 
locations  of  the  radius  of  the  central  jet,  inner  and  outer  radii 
of  the  swirling  stream,  and  the  radius  of  the  centerbody  are 
identified  by  diamond  markers  shown  along  the  ordinate. 


Profiles 


Turbulent  Klneti c  Enerc 


Figures  8a  through  8h  show  the  radial  variation  of 
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(dimensional)  turbulent  kinetic  energy  [expressed  in  (m/s)  ]  for 
the  eight  test  cases.  Note  that  with  the  assumption  of  isotropy, 
these  profiles  also  represent  the  variation  of  the  square  of  the 
rms  axial  velocity  components  (magnified  by  a  factor  of  3/2). 


These  profiles  show  that  in  all  the  cases  relatively 
strong  local  peaks  occur  in  the  vicinity  of  the  central-jet 
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radius.  This  signifies  the  development  of  the  central  Jet. 
Previous  LDA  measurements  (Reference  10)  of  the  development  of 
the  central  Jet  in  a  weak  annular  flow  also  observed  similar 
trends.  The  local  peaks  in  the  shear  layers  of  the  inner  and 
outer  streams,  on  the  other  hand,  do  not  exhibit  a  relative 
dominance,  unless  the  flowfields  are  characterized  by  increased 
strength  of  these  streams.  Thus,  Figure  8c  (corresponding  to 
case  3)  shows  the  greatly  enhanced  turbulence  activity  in  the 
shear  layers  of  the  swirling  stream.  As  discussed  in  Paragraph 
Ill.l.a(l),  case  3  is  an  example  of  the  flowfield  dominated  by 
the  swirling  stream  (recall  that  in  case  3,  ms  is  twice  that  of 
case  1  and  the  swirl  number  is  the  highest) .  An  integrated 
perspective  of  the  increased  turbulent  mixing  in  the  near  wake  is 
available  by  inspecting  Figures  3c  and  8c. 

The  relative  dominance  by  the  outer  annular  stream 
is  noticeable  in  Figures  8g  and  8h  (which  correspond  to  cases  7 
and  8  with  m  =  2  kg/s)  which  show  the  increasing  peak  of 
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turbulent  kinetic  energy  in  the  outer  shear  layer  (at  approxi¬ 
mately  one  radius  of  the  centerbody) .  Finally,  the  effect  of 
swirl  in  the  flowfield  can  be  understood  by  comparing  Figures  8f 
and  8h  with  Figures  8e  and  8g  (swirl-free  cases)  respectively. 

The  introduction  of  swirl  serves  to  increase  the  turbulence 
activity  in  the  vicinity  of  the  middle  stream. 

It  is  of  interest  to  compare  these  results  with 
respect  to  the  predictions  in  the  old  centerbody  configuration 
(Reference  3).  The  radial  distribution  of  the  turbulent  kinetic 
energy  for  2  kg/s  annular  flow  and  zero  central  flow  (seen  in 
Figure  29  of  Reference  3)  shows  trends  similar  to  those  in  Figure 
8g .  The  comparison,  however,  of  the  predicted  and  measured 
radial  distributions  of  the  rms  axial  velocity  field  in  a 
small-scale  combustor  (seen  in  Figure  19  of  Reference  3)  suggests 
that  in  the  case  of  present  predictions  also,  the  quantitative 
comparison  with  measurements  (when  they  are  available)  is  likely 
to  be  poor  [see  Paragraph  III . 1 . b( 3b) ] . 


Figures  9a  through  9h  show  the  radial  variation  of 
the  (dimensional)  temperature  field  at  four  social  stations.  With 
our  choice  of  the  initial  temperatures  of  400°K  for  the  central 
jet  and  293 °K  for  the  two  annular  streams  and  with  the  rather 
small  mass-flow  rate  m  ,  it  should  not  be  surprising  that  the 
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temperature  profiles  indicate  only  a  very  small  spatial  extent  of 
the  temperature  rise  above  the  ambient.  Such  dramatic  effects  as 
seen  for  the  velocity  fields  (e.g.,  cases  3,  7,  and  8)  are  not 
discernible  for  the  temperature  field.  Perhaps  much  hotter 
central  jet,  or  an  increase  in  the  temperature  of  the  swirling 
stream  could  show  more  interesting  characteristics  of  the 
temperature  field  as  it  is  affected  by  the  different  parameters 
of  the  problem. 


( 3 )  Profiles  of  C02  Mass  Fraction  and  Swirl  Velocity 

The  discussion  of  the  radial  distributions  is 
concluded  with  the  consideration  here  of  the  profiles  of  the 
C02  mass  fraction  and  the  mean  swirl  velocity  fields.  These  are 
seen  in  Figures  10a  through  lOf.  The  four  axial  stations 
considered  here  represent  locations  much  farther  downstream  (by  a 
factor  of  5  to  10)  than  those  seen  earlier  for  the  profiles  of 
temperature  and  kinetic  energy.  The  present  choice  is  dictated 
more  by  the  increased  penetration  of  the  flowfield  by  the  mean 
swirl  velocity  component  in  both  axial  and  radial  directions. 

In  the  figures  the  C02  mass  fraction  profiles  correspond  to  those 
that  are  essentially  confined  to  small  values  of  the  ordinate. 
Note,  however,  that  the  C02  mass  fraction  profiles  will  be 
identical  to  the  temperature  [normalized  as  in  Paragraph 
III.l.b(2)]  profiles  seen  in  Figures  9a  through  9h  (as  observed 
previously  in  Figures  6a  through  6h) .  Therefore,  attention  in  the 
present  discussion  is  essentially  directed  to  an  examination  of 
the  swirl  velocity  profiles. 
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Figure  9(b).  Radial  Profiles  of  Temperature. 
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Figure  10(c).  Redial  Profile*  of  CO  Mas*  Fraction  and  Swirl 
Velocity. 
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Figure  10(d) .  Radial  Profiles  of  CO.  Hass  Fraction  and  Swirl 
Velocity.  * 
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The  distributions  of  the  swirl  velocity  in  the 
figures  correspond  to  those  profiles  which  display  the 
of f centerline  (R/H  =  0)  peaks.  The  peak  values  of  the  swirl 
velocity  occur  in  the  vicinity  of  the  swirling  stream  at  the  two 
axial  stations  closest  to  the  bluff  body.  This  is  indicative  of 
the  fact  that  the  swirl  remains  confined  to  the  middle  stream. 

The  peaks  are  shifted  radially  outward,  towards  the  outer  annular 
stream,  at  the  two  stations  farther  downstream.  It  is  interest¬ 
ing  to  note  that  all  these  peaks  occur  on  the  outer  side  of  the 
swirling  stream  in  all  the  cases.  That  no  peak  occurred  on  the 
inner  side  of  the  swirling  stream  in  all  the  parametric  vari¬ 
ations  considered  by  us  suggests  that  this  is  more  likely  to  be 
a  consequence  of  the  present  geometrical  parameters.  This  con¬ 
jecture  can  be  ascertained  only  by  a  CPD  examination  of  a 
parametric  variation  of  the  swirler  and  centerbody  radii. 

The  profiles  at  the  two  near-field  locations 
conform  to  each  other  more  closely  (this  is  especially  true  in 
Figure  lOe) .  The  more  dramatic  effect  of  the  stronger  swirling 
stream  is  brought  out  in  Figure  10c  which  shows  the  relative 
increase  in  swirl  velocity  even  at  the  two  far-field  locations. 
Another  interesting  aspect  of  the  profiles  at  these  two  latter 
locations  is  noticeable  in  Figures  10c,  lOd,  and  lOe.  The  radial 
profiles  in  these  cases  show  a  much  greater  symmetry,  unlike  the 
clearly  asymmetric  profiles  (more  skewed  towards  the  outer 
annular  stream)  seen  typically.  The  physical  significance  of 
this  peculiar  behavior  is  not  clear.  We  must  determine  from 
experiments  how  realistic  is  the  predicted  behavior  concerning 
this  and  other  aspects  of  the  flowfield. 

Although  the  radial  distributions  discussed  herein 
have  not  addressed  the  mean  axial  velocity  fields,  it  is 
generally  likely  that  these  profiles  will  exhibit  trends  observed 
in  the  previous  POSF  configuration.  The  agreement  with 
experimental  data  is  also  likely  to  be  fair.  Finally,  we  may 
expect  that  Abramovich-type  universal  profiles  (see  Reference  2) 
will  also  apply  in  the  proposed  configuration. 


2.  REACTING  FLOWFIELDS 


Cases  9  and  10  of  the  present  CFD  study  are  concerned  with 
the  reacting  flowfields  of  the  proposed  POSF  configuration  which 
result  from  the  combustion  of  gaseous  propane  in  the  central  jet 
and  air  in  the  two  annular  streams.  Of  the  two  cases,  case  9 
deals  with  the  swirl-free  middle  stream,  while  case  10 
corresponds  to  the  swirling  middle  stream.  Before  we  discuss 
each  case  in  turn,  some  comments  are  in  order  regarding  the 
reacting-f lowf ield  calculation  procedure  of  Reference  5. 

a.  Shortcomings  of  Computational  Procedure 

The  computer  code  of  Reference  5  employs  the 

Magnussen-H jertager  model  (Reference  7)  for  describing  the 

combustion  of  propane  and  air  (see  Paragraph  II.  2).  The  source 

term  S  for  the  fuel-mass-fraction  equation  [see  Equation  (1)], 
mf 

in  accordance  with  this  model,  is  explicitly  indicated  in  Table 
1 .  The  appropriate  value  of  Smf  is  denoted  by  the  minimum  of  the 
three  terms  found  inside  the  square  brackets.  While  a  strict 
implementation  of  the  Magnussen-Hjertager  model  would  require  the 
search  for  this  minimum  value,  the  computer  code  of  Reference  5 
has  not  complied  with  that  requirement.  The  reason  for  this 
lacuna  in  the  code  is  not  clear.  When  the  present  study  carried 
out  a  strict  Implementation  of  the  Magnussen-Hjertager  model  In 
the  calculations,  the  computed  results  led  to  nonphysical  and 
large  negative  mass  fractions  of  oxygen.  Presumably,  this 
difficulty  might  have  persuaded  the  code  developers  to  abandon 
the  correct  implementation.  For  the  present,  it  is  essential  to 
recognize  this  deficiency  of  the  computational  procedure. 


Furthermore,  the  computer  code  Imposes  an  arbitrary  lower 
bound  on  the  temperature  field.  The  prescribed  minimum  value  is 
293*K  (see  line  83  of  subroutine  PROPS  of  the  code).  When  this 
was  changed  to  150*K,  we  found  that  the  temperatures  in  the 
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flowfield  reached  values  as  low  as  160#K.  It  appears  that  major 
modifications  are  required  in  the  algorithm  to  implement  a  proper 
coupling  of  the  combustion  model  to  the  hydrodynamics.  Clearly, 
the  required  algorithm  improvements  are  beyond  the  scope  of  the 
present  research  study.  Instead,  the  present  intent  has  been  to 
examine  two  computational  cases,  mainly  to  obtain  some  tentative 
ideas  on  the  reacting  flowfield  behavior  of  the  new  POSF 
configuration . 

b .  Swirl-Free  Flowfield 

Case  9  has  a  value  of  zero  for  the  parameter  S_ . 

The  mass-flow  rates  in  all  the  three  streams  are  identical  with 
the  corresponding  values  of  cases  7  and  8  of  the  nonreacting 
flows.  Thus,  this  flowfield  is  largely  an  annular-stream 
dominant  one.  Note  that  a  comparison  of  cases  7  and  9  will 
reveal  the  influence  of  combustion.  Before  we  look  at  the 
computed  results,  it  is  important  to  mention  that  the  present 
calculations  never  converged  for  case  9.  While  convergence  of 
the  numerical  calculations  depends  on  the  specific  value  of  the 
many  inflow  parameters  of  the  problem,  no  attempt  was  made  to 
investigate  the  parametric  effect  on  the  numerical  convergence 
for  this  particular  case. 

( 1 )  Overview  of  Flowfield  Structure 

Figures  11a  through  11c  show  the  plots  of 
velocity  vectors,  temperature  contours,  and  density  contours 
respectively.  An  overview  of  the  flow  structure  can  be  obtained 
from  an  inspection  of  these  three  plots.  Furthermore,  the 
changes  in  the  flowfield  caused  by  the  introduction  of  combustion 
are  noted  by  comparing  Figure  11a  with  Figure  4b  (the  bottom 
half) .  Thus,  combustion  is  seen  to  generate  an  additional  vortex 
in  Figure  11a.  The  vortex  center  of  this  clockwise  vortex  is 
located  farther  downstream  of  the  primary  vortex  center  and 
closer  to  the  centerline  than  the  latter.  It  appears  as  if  the 


Figure  11(b).  Temperature  Contours  for  Swirl-Free  Reacting  Flow 


Figure  11(c).  Density  Contours  for  Swirl-Free  Reacting  Flow. 


primary  vortex  in  the  nonreacting  case  (Figure  4b)  has  split  into 
two  large  vortices  moving  away  from  each  other.  Their  spatial 
separation  is  large  enough  that  no  counterclockwise  vortex  is 
required  in  between  (on  kinematic  grounds) .  Their  separation  is 
enhanced  by  a  further  penetration  of  the  central  jet.  This  is 
evident  from  the  downstream  movement  of  the  forward  stagnation 
point  to  a  normalized  distance  (in  terms  of  the  centerbody 
diameter)  of  0.83  which  is  nearly  twice  the  nonreacting  value  of 
0.465.  The  downward  movement  of  the  rear  stagnation  point  is 
slight  however  (from  1.455  previously  to  1.475  now).  The  marked 
increase  in  the  forward  stagnation  point  location  due  to 
combustion  is  consistent  with  the  experience  in  the  old 
centerbody  configuration  (Reference  19) . 

Further  elucidation  of  the  reacting  flowfield  can  be 
had  from  Figures  lib  and  11c.  The  regions  of  the  maximum 
temperature  and  minimum  density  in  the  flowfield  can  be 
identified  by  the  dark  bands  in  the  temperature  and  density 
contours.  Clearly  these  regions  correspond  to  the  locations 
where  the  combustion  intensity  of  the  flame  is  a  maximum.  A  look 
at  Figure  11a  shows  that  the  flame  occurs  in  the  region  where  the 
fuel  is  drawn  in  the  counterclockwise  direction  and  is  mixed  with 
the  incoming  fresh  air  of  the  middle  stream  first  and  the 
entrained  air  of  the  outer  annular  stream  (through  the  clockwise 
vortex)  next.  The  supply  of  fuel  into  the  reaction  zone  is 
augmented  both  by  an  entrainment  process  due  to  a  small 
counterclockwise  vortex  (located  between  the  central  jet  and  the 
middle  stream)  and  by  the  flow  reversal  and  forward  stagnation  of 
the  central  jet  due  to  the  bigger  clockwise  vortex  present 
farther  downstream.  Note  that  the  latter  also  aids  in  bringing 
the  hot  combustion  products  into  the  reaction  zone,  thereby 
providing  the  required  species  and  energy  for  the  chemical 
reactions . 


{ 2 )  Centerline  Variations 

The  foregoing  description  of  the  overall  flowfield 
facilitates  the  considerations  of  the  centerline  variations  of 
the  axial  velocity  fields  and  the  scalar  fields.  These  are 
presented  in  Figures  12a  and  12b  respectively.  The  normalized 
mean  axial  velocity  profile  in  Figure  12a  may  be  compared  to  the 
corresponding  nonreacting  case  seen  in  Figure  5g.  Such  a 
comparison  immediately  reveals  that  with  the  occurrence  of 
combustion,  the  mean  axial  velocity  decays  more  slowly;  the 
central  jet  penetrates  farther;  the  forward  stagnation  point  is 
pushed  farther  downstream  (nearly  twice  as  far  as  in  Figure  5g) ; 
the  magnitude  of  the  peak  negative  velocity  is  decreased  by  more 
than  a  factor  of  2  and  its  location  is  moved  farther  downstream; 
and  the  location  of  the  rear  stagnation  point  and  the  subsequent 
recovery  of  the  mean  axial  velocity  are  only  slightly  affected. 
All  these  trends  are  generally  consistent  with  the  experience  in 
the  earlier  POSF  configuration  (e.g.,  see  Reference  19). 

Comparison  of  the  nonreacting  and  reacting  profiles 
of  the  rms  velocity  component  shows  that  there  are  slight 
differences  in  the  magnitudes  and  locations  of  the  two  peaks. 
Overall,  however,  the  turbulence  intensity  in  the  reacting  case 
is  seen  to  be  higher  almost  everywhere  along  the  centerline. 
While  the  assumption  of  isotropy  has  a  crucial  effect  on  this 
prediction,  it  appears  that  combustion  in  the  present  situation 
of  initially  unmixed  reactants  might  offset  the  decrease  in 
turbulence  due  to  dilatation  effect  by  an  increase  due  to  the 
production  of  turbulence  kinetic  energy  through  augmented  shear 
stress.  Future  measurements  should  be  able  to  verify  this 
possibility. 

The  scalar  profiles  in  Figure  12b  show  the 
centerline  variations  of  the  mass  fraction  of  propane  and  the 
normalized  temperature  (T  -  Tmin)/(Tmax  “  Tmin) *  Here  Tmax  is 
the  flame  temperature.  The  scalar  profiles  show  that  the 
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centerline  temperature  peaks  where  the  propane  mass  fraction 
vanishes.  Furthermore,  from  Figure  12a  this  location  is  seen  to 
coincide  with  the  forward  stagnation  point.  Note,  however,  that 
the  centerline  peak  temperature  is  only  83*  of  the  normalized 
flame  temperature.  This  confirms  that  the  maximum  temperature 
(i.e.,  the  most  intense  reaction  zone)  occurs  of f centerline.  We 
note  that  the  flame  temperature  in  case  9  is  1,850°K. 

( 3 )  Radial  Profiles 

Before  concluding  the  discussion  of  case  9,  we 
consider  the  radial  variation  of  turbulent  kinetic  energy  and 
temperature  at  different  axial  locations.  These  distributions 
(for  dimensional  values  of  the  variables)  are  presented  in 
Figures  13a  and  13b  respectively. 

A  comparison  of  Figures  8g  and  13a  shows  that  there 
is  very  little  difference  in  the  trends  of  the  radial  profiles  of 
turbulence  kinetic  energy  between  the  nonreacting  and  reacting 
cases.  Indeed,  in  the  two  near-field  locations  the  profiles 
almost  coincide.  In  the  two  far-field  locations,  the  profiles 
coincide  in  the  outer  annular  region,  indicating  thereby  the 
very  little  influence  of  combustion  on  the  turbulence  there. 

Closer  to  the  centerline,  the  kinetic  energy  In  the  reacting  case 
is  more  than  that  in  the  nonreacting  case,  the  increase  being 
appreciable  especially  at  the  farthest  axial  location.  These 
flowfield  features  are  in  conformity  with  the  aspects  discussed 
earlier . 

The  temperature  profiles  seen  in  Figure  13b  confirm 
the  overall  description  suggested  by  the  velocity-vector  plot  and 
temperature-contour  plot  seen  in  Figures  11a  and  lib  respectively. 
The  propagation  of  the  reaction  zone  axially  downstream  and 
radially  outward  is  clearly  noticeable  from  the  radial 
distributions  of  the  temperature.  Also  unmistakable  is  the 
negligible  impact  of  combustion  on  the  outer  parts  of  the 
f lowf ield . 


TURBULENT  KINETIC  ENERGY 


ire  13(b).  Radial  Profiles  of  Temperature 
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c .  Flowfield  with  Swirl 

The  flowfield  with  nonzero  swirl  in  the  middle  stream 
(SR  =  1*0)  corresponds  to  case  10.  Note  that  the  computations  in 
this  case,  unlike  those  of  case  9,  had  converged.  The  precise 
reason  for  the  favorable  influence  of  the  presence  of  swirl  on 
the  numerical  convergence  is  not  clear. 

( 1 )  Overall  Flowfield  Features 

The  global  features  of  the  flowfield  can  be 
understood  from  the  plots  of  velocity  vectors,  temperature 
contours,  and  density  contours  seen  in  Figures  14a,  14b,  and  14c 
respectively.  Furthermore,  a  comparison  of  cases  9  and  10  shows 
the  effect  of  swirl  on  the  reacting  flowfield,  whereas  a 
comparison  of  cases  8  and  10  reveals  the  influence  of  combustion 
on  the  swirling  flowfield. 

From  a  comparison  of  Figures  11a  and  14a  it  is  clear 
that  while  the  overall  flowfield  structure  remains  the  same  with 
or  without  swirl,  the  details  concerning  the  location  of  the 
vortex  centers  and  the  spatial  extent  of  the  vortices  are 
different.  Thus,  with  swirl  the  large  clockwise  vortex  farthest 
from  the  bluff  body  has  moved  radially  outward  and  axially 
upstream.  This  movement  towards  the  centerbody  shoulder  in  turn 
causes  the  upstream  movement  of  the  other  large  clockwise  vortex, 
as  well  as  the  small  counterclockwise  vortex  (the  spatial  extent 
of  the  latter  is  almost  halved  with  swirl) .  These  trends  were 
also  seen  [in  Paragraph  III.l.a(2)]  for  nonreacting  flows  (cases 
7  and  8).  The  small  counterclockwise  vortex  (located  between  the 
swirling  stream  and  the  central  Jet),  on  the  other  hand,  is  seen 
to  have  moved  axially  downstream  and  radially  outward  (with  the 
radial  movement  being  much  smaller  than  the  axial  one)  when  swirl 
is  present.  Qualitatively  speaking,  the  net  result  of  these 
vortex  movements  is  to  intensify  the  mixing  rate  in  the  L-shaped 
region  formed  by  the  three  leftmost  vortices.  This  is  also 
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Figure  14(b).  Temperature  Contours  for  Reacting  Flow  with  Swirl. 


Figure  14(c).  Density  Contours  for  Reacting  Flow  with  Swirl. 


evident  from  the  more  tightly  bound  system  of  vortices  in  Figure 
14a,  in  comparison  with  that  in  Figure  11a.  More  quantitatively, 
the  normalized  locations  of  the  forward  and  rear  stagnation 
points  with  swirl  are  0*735  and  1*62  (the  former  has  moved  up¬ 
stream  and  the  latter  downstream) .  The  enhancement  of  the  mixing 
and  combustion  processes  resulting  from  the  improved  flowfield 
characteristics  due  to  the  introduction  of  swirl  is  demonstrated 
by  the  flame  temperature  increasing  to  1,905°K  (from  1850*K  for 
case  9) . 

The  changes  in  the  swirling  flowfield  due  to  the 
introduction  of  combustion  [seen  from  a  comparison  of  Figures  14a 
and  4b  (the  top  half)  corresponding  to  case  8  and  10]  are 
similar  to  the  changes  noted  in  the  nonswirling  flowfield  due  to 
combustion  [see  Paragraph  III.2.b(l)  for  the  comparison  of  cases 
7  and  9] .  Thus,  we  encounter  all  the  features  such  as  the 
splitting  of  the  large  clockwise  outer  vortex  into  two  clockwise 
vortices  (that  are  spatially  separated) ,  the  greater  penetration 
of  the  central  jet,  and  the  resulting  increase  in  the  normalized 
distances  of  the  centerline  forward  and  rear  stagnation  points 
(the  latter  to  a  lesser  extent).  A  quantitative  comparison  can 
be  easily  established  from  Figures  14a  and  4b  (the  top  half). 

Further  insights  on  the  reacting  flowfield  are 
available  from  the  contours  of  the  scalar  fields  in  Figures  14b 
and  c.  A  comparison  of  Figures  lib  and  14b  shows  that  the 
introduction  of  swirl  serves  to  confine  the  reaction  zone  closer 
to  the  bluff  body  and  that  the  overall  effect  is  to  enhance 
combustion. 

( 2 )  Centerline  Variations 

The  centerline  variations  of  the  axial  velocity 
components  (mean  and  rms),  the  scalars  (propane  mass  fraction  and 
temperature) ,  and  the  mean  swirl  velocity  component  are  shown  in 
Figures  15a-c  respectively.  The  effect  of  swirl  on  the  reacting 
flow  for  the  axial  velocity  fields  and  the  scalar  fields  can  be 
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obtained  from  a  respective  comparison  of  Figures  12a  and  15a,  and 
12b  and  15b.  The  changes  in  both  the  mean  and  rms  components  are 
not  appreciable.  There  are  slight  increases  in  both  stagnation 
point  distances,  as  well  as  in  the  magnitude  and  location  of  the 
peak  negative  mean  velocity.  There  is  a  small  decrease  in  the 
first  peak  and  a  small  increase  in  the  second  peak  in  the  rms 
profiles . 

These  observations  are  confirmed  for  the  scalar 
profiles  also.  There  is  a  very  little  change  in  the  propane  mass 
fraction  profiles  in  Figures  12b  and  15b.  The  temperature 
profiles  exhibit  the  same  overall  trends.  However,  the  peak 
centerline  value  of  the  normalized  temperature  in  Figure  15b  is 
only  77*.  This  shows  that  while  the  flame  temperature  has 
increased  by  about  55®  with  the  introduction  of  swirl,  the  peak 
centerline  value  has  dropped  by  51°  [0*77  (1905-293)  +  293  vs. 
0*83  (1850-293)  +  293].  This  appears  to  be  consistent  with  the 
greater  confinement  and  intensification  of  mixing  and  combustion 
in  the  near-wake  of f centerline  regions  caused  by  the  introduction 
of  swirl. 

The  effect  of  combustion  on  the  centerline  profiles 
of  the  swirling  flowfield  can  be  easily  understood  by  a  compari¬ 
son  of  Figures  5h  and  15a,  6h  and  15b,  and  7f  and  15c  for  the 
axial  velocity  components,  scalar  fields,  and  the  mean  swirl 
velocity  respectively. 

As  seen  earlier  for  the  swirl-free  flowfield  [in 
Paragraph  III.2.b(2)],  the  effect  of  combustion  on  the  mean  and 
rms  velocity  components  is  to  produce  similar  trends,  viz., 
greater  penetration  by  the  central  jet,  increased  centerline 
stagnation  point  distances,  slower  decay  of  the  mean  velocity, 
and  an  overall  increase  almost  everywhere  of  the  rms  velocity. 
Thus,  all  our  earlier  comments  in  connection  with  the  nonswirling 
flow  apply  here  also.  An  interesting  exception,  however,  is  that 
while  the  overall  rms  velocity  trends  are  preserved  in  cases  7 


and  9  (Figures  5g  and  12a),  In  case  10  the  effect  of  combustion 
is  much  more  pronounced  on  the  first  rms  velocity  peak.  Thus, 
whereas  the  first  peak  Is  smaller  than  the  second  in  Figure  5h 
(case  8),  it  becomes  larger  in  Figure  15a  (case  10).  This 
behavior  (which  will  be  further  clarified  in  our  discussion  of 
the  radial  profiles)  provides  additional  evidence  of  the  greater 
turbulence  activity  in  the  near  field. 

In  the  comparison  of  the  scalar  fields,  very  little 
of  significance  could  be  said,  in  view  of  the  direct  and  strong 
effect  of  combustion  on  both  the  temperature  and  propane  mass 
fraction.  Note,  however,  that  the  Increased  penetration  by  the 
central  jet  in  the  reacting  flowfleld  is  clearly  evident  from  the 
slower  decay  of  the  propane  mass  fraction  (in  Figure  15b)  when 
compared  to  the  decay  of  CO 2  mass  fraction  (in  Figure  6h) .  Of 
course,  unlike  case  8  where  the  profiles  of  C02  mass  fraction  and 
temperature  coincide,  case  10  properly  describes  the  peaking  of 
the  temperature  caused  by  the  exothermlcity  of  combustion. 

Figure  15c  and  7f  bring  out  the  effect  of  combustion 
on  the  mean  swirl  velocity  variation  close  to  the  centerline. 
Recall  the  earlier  discussion  [in  Paragraph  III.l.b(4)],  however, 
regarding  the  negligibly  small  magnitudes  of  the  swirl  velocity 
in  the  vicinity  of  the  centerline.  Keeping  this  in  mind,  we  see 
that  the  normalized  swirl  velocity  profile  differs  negligibly 
initially.  Although  the  peak  value  has  Increased  by  nearly  30% 
with  combustion,  the  location  of  this  peak  has  only  slightly 
shifted  downstream.  A  more  significant  change  is  noted 
subsequently,  as  the  proximity  to  the  reaction  zone  Increases. 

The  swirl  velocity  begins  to  decay  much  more  slowly  and  when  the 
local  minimum  occurs,  its  location  is  approximately  twice  as  far 
(in  Figure  15c)  as  that  in  the  nonreacting  case  (Figure  7f ) . 

This  behavior  is  consistent  with  the  Increased  penetration  by  the 
central  jet.  Farther  downstream  of  the  minimum  location,  there 
is  no  appreciable  effect  of  combustion  on  the  swirl  velocity. 
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(3)  Radial  Profiles 

Finally,  we  examine  the  radial  distributions  of 
turbulence  kinetic  energy,  temperature,  propane  mass  fraction, 
and  mean  swirl  velocity  before  concluding  the  discussion  of  case 
10.  Figures  16a  through  16c  present  these  radial  profiles. 

The  swirl  in  the  middle  stream  does  not  appear  to 
have  a  significant  effect  on  the  turbulence  kinetic  energy 
profiles.  This  is  readily  apparent  by  comparing  Figure  16a  with 
Figure  13a.  Only  at  the  farthest  axial  station  is  there  some 
redistribution  of  the  turbulence  energy.  With  swirl,  the 
magnitudes  of  the  peaks  at  the  extremities  (i.e.,  near  the 
centerline  and  towards  the  outer  annular  stream)  are  slightly 
decreased  and  the  magnitude  in  the  middle  (from  the  inner  radius 
of  the  swirling  stream  to  about  eight-tenths  of  centerbody 
radius)  is  increased.  At  the  two  axial  stations  farther  upstream 
too,  the  incipient  tendency  of  swirl  to  increase  the  turbulence 
kinetic  energy  in  the  middle  can  be  noticed,  but  the  increase  is 
not  significant.  Overall,  it  can  be  safely  said  that  the  intro¬ 
duction  of  swirl  certainly  causes  an  increase  in  the  turbulence 
activity  in  the  mid-radial  regions,  an  outcome  that  contributes 
to,  and  is  consistent  with,  the  enhanced  mixing  and  combustion 
alluded  to  earlier. 

The  effect  of  combustion  on  turbulence  kinetic 
energy  can  be  discerned  from  comparing  Figure  16a  with  Figure  8h. 
Minor  differences  are  noted  between  the  two  cases,  especially  at 
the  peaks.  Below  the  inner  radius  of  the  swirling  stream,  some 
increase  in  the  kinetic  energy  is  seen  at  the  second  and  third 
axial  stations.  At  the  farthest  station,  an  increase  in  the 
near-centerline  peak  and  a  decrease  in  the  outer  peak  are  also 
noted.  These  trends  are  similar  to  those  discussed  for  non- 
swirling  flows  in  Paragraph  III.2.b(3)  and  are  consistent  with 
our  anticipations. 
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TURBULENT  KINETIC  ENERGY 
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Figure  16b  shows  the  radial  profiles  of  the 
dimensional  temperature  at  the  four  axial  stations.  As  seen  In 
Figure  13b  for  the  swlrl-free  reacting  flow,  the  profiles  in 
Figure  16b  also  display  the  propagation  of  the  reaction  zone  in  a 
radially  outward  and  axially  downstream  direction.  A  respective 
comparison  of  the  individual  plots  In  Figures  13b  and  16b  reveals 
the  effect  of  swirl  on  the  radial  temperature  profiles.  At  the 
two  near-field  locations,  the  profiles  show  that  the  only 
difference  in  them  is  a  decrease  in  the  peaks  by  about  20ft  due 
to  the  introduction  of  swirl.  Farther  downstream,  this  decrease 
begins  to  get  smaller  and  there  is  also  a  radial  spread  of  the 
hotter  regions.  This  is  dramatically  established  at  the  farthest 
station  where  the  peak  temperature  is  larger  with  swirl;  the 
radial  spread  of  the  hotter  regions  is  nearly  doubled;  and  the 
region  between  the  inner  radius  of  the  swirling  stream  and  eight- 
tenths  of  centerbody  radius  is  characterized  by  uniformly  high 
temperatures  (exceeding  1600*K) .  This  particular  profile  clearly 
demonstrates  that  the  introduction  of  swirl  greatly  enhances  the 
combustion  and  heat  release  processes  in  the  mid-radial  regions. 

The  radial  distributions  of  propane  mass  fraction 
and  mean  swirl  velocity  are  shown  in  Figure  16c.  As  seen  in 
Figure  10  (a  through  f)  of  nonreacting  flowfields,  the  four 
profiles  of  the  swirl  velocity  are  those  which  exhibit  the 
of f centerline  peaks.  The  propane  mass  fraction  profiles  are 
monotonlc  and  closer  to  R/H  *0.  At  the  two  near-field  stations, 
the  radial  spread  of  propane  is  small  and  there  is  negligible 
mixing  of  air  at  these  radial  locations.  At  the  third  axial 
station,  nonzero  values  of  propane  mass  fraction  can  be  seen  as 
far  as  eight-tenths  of  centerbody  radius.  Note  too  the  signi¬ 
ficant  mixing  with  air  in  these  radial  regions  (e.g.,  the  maximum 
value  of  propane  mass  fraction  is  only  38ft  at  R/H  «  0).  At  the 
farthest  axial  station,  of  course,  propane  is  no  longer  present, 
having  been  consumed  farther  upstream  at  all  radial  locations. 

An  inspection  of  Figure  iOf  (for  the  nonreacting  case  8)  shows 
that  the  radial  spread  of  C02  is  much  less.  Also  at  the  two 


far-field  stations,  nonzero  values  of  C02  mass  fraction  were 
still  present.  Therefore,  the  introduction  of  swirl  is  seen  to 
have  improved  the  mixing  of  propane  and  air  in  the  near-field 
of f centerline  locations  and  thus  to  have  contributed  to  better 
and  more  complete  combustion.  This  observation,  of  course,  is 
consistent  with  the  radial  temperature  profiles  discussed 
previously  (note  the  different  axial  locations  in  Figures  16b  and 
16c)  . 

Finally,  the  swirl  velocity  profiles  in  Figure  16c 
show  the  effect  of  combustion  when  they  are  compared  with  the 
corresponding  profiles  of  Figure  lOf.  We  note  that  at  the  two 
near-field  axial  locations  there  are  minor  differences.  The  peak 
normalized  values  at  the  nearest  location  are  almost  unchanged  in 
their  magnitudes  (note  the  difference  in  the  scales)  and  radial 
locations.  At  other  radial  locations  above  and  below  the  peaks, 
the  swirl  velocity  with  combustion  is  generally  higher.  This  is 
especially  true  of  the  smaller  peaks  closer  to  the  centerline. 
Although  the  radial  (outward)  movement  of  this  peak  is  small  with 
combustion,  its  magnitude  is  almost  trebled.  This  behavior  of 
enhanced  swirl  velocity  due  to  combustion  is  consistent  with  the 
earlier  observations.  These  features  are  also  exhibited  by  the 
profile  at  the  second  near-field  station.  In  addition,  while  the 
outward  spread  of  swirl  remains  essentially  the  same,  the  inward 
spread  (i.e.,  towards  the  centerline)  of  swirl  with  combustion 
has  increased.  No  doubt,  this  has  contributed  to  the  improved 
combustion  in  these  radial  regions.  At  the  two  far-field 
locations,  the  magnitudes  of  the  swirl  velocity  are  quite  small. 
But  even  in  these  locations,  the  Increase  due  to  combustion  is 
noticeable.  For  Instance,  the  swirl  velocity  at  eight-tenths  of 
centerbody  radius  is  doubled  with  combustion. 


SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  section  outlines  the  major  conclusions  emerging  from  the 
CFD  investigations  of  the  present  Scholarly  Research  Program  and, 
in  the  light  of  these  conclusions,  offers  some  recommendations 
concerning  further  assessment  of  computations  of  gas  turbine 
combustor- like  flowfields. 

1 .  CONCLUSIONS 

Numerical  predictions  of  the  recirculating  confined  turbulent 
flowfields  with  and  without  swirl  in  both  nonreacting  and  react¬ 
ing  situations  have  been  made  for  the  proposed  POSF  combustor 
configuration  under  the  framework  of  t ime- independent ,  Reynolds- 
averaged  Navier-Stokes  equations.  No  time-dependent  formulation 
has  been  considered  within  the  scope  of  the  present  program  for 
the  prediction  of  the  proposed  configuration,  either  for 
obtaining  the  steady-state  solutions,  or  for  describing  the 
dynamic  flowfield  features  likely  to  be  encountered.  However,  an 
assessment  of  the  solutions  of  the  time-dependent,  Reynolds- 
averaged  Navier-Stokes  equations  for  the  nonreacting  flowfields 
in  the  existing  POSF  centerbody  configuration  through  the 
MacCormack  algorithm  is  documented  in  Volume  I  of  this  report. 
Some  of  the  conclusions  and  recommendations  reported  there  must 
be  considered  in  conjunction  with  those  documented  here. 

Following  are  the  main  conclusions  of  this  study: 

•  Calculations  using  the  t ime- Independent ,  Reynolds-aver¬ 
aged  formulation  with  the  k-e  turbulence  model  have 
furnished  the  numerical  predictions  of  the  complex 
turbulent  flowfields  in  the  proposed  ducted  bluff-body 
combustor  configuration. 


A  limited  parametric  examination  of  the  various  fluid 
dynamical  aspects  relating  to  the  two  annular  air  streams 


(swirl-free  outer  and  swirling  inner  flows)  and  the 
central  fuel  jet  has  addressed  the  effects  of  aass  flow 
rates,  swirl,  and  combustion  on  the  computed  flowfields. 

The  CPD  examination  of  the  steady-state  flowfields  has 
considered  the  features  of  the  global  flowfleld  from  the 
velocity-vector  plots  (and  also  from  the  contours  of  the 
temperature  and  density  fields  in  reacting  flows),  the 
variations  of  the  centerline  profiles  of  selected 
variables,  and  the  radial  distributions  of  these 
variables  at  a  few  axial  locations;  and  thereby  obtained 
descriptions  of  the  complex  flowfleld  Interactions 
Inherent  in  the  proposed  configuration. 

The  success  of  the  present  CFD  predictions  remains  to  be 
verified  by  experimental  measurements  in  the  POSF 
combustor.  Thus,  the  present  testing  of  CFD  methods 
against  hitherto  nonexistent  and  nonmeasured  flowfields 
represents  true  predictions  whose  validity  is  ascertain¬ 
able  by  experimental  means. 

Although  the  present  predictions  remain  to  be  validated, 
the  Internal  consistency  exhibited  by  these  predictions 
and  their  conformity  with  some  of  the  predicted  behavior 
verified  in  the  existing  centerbody  configuration  encour¬ 
age  the  belief  that  the  numerical  predictions  will  show 
reasonable  agreement  with  the  actual  flowfleld  behavior. 

The  success  anticipated  for  the  predictions  in  conforming 
to  the  realistic  behavior,  however,  appears  to  be  re¬ 
stricted  to  the  mean-field  description.  The  quantitative 
agreement  of  the  predictions  of  the  rms  fields  with  the 
experimental  data  is  likely  to  remain  poor.  Again,  even 
in  the  mean  field,  the  anticipated  success  in  the  react¬ 
ing-flow  predictions  is  unlikely  to  match  that  for  non- 
reacting-  flow  predictions. 


•  Nevertheless,  the  present  predictions  are  expected  to 
offer  significant  insights  into  the  development  and 
testing  of  the  proposed  POSF  combustor  to  simulate  some 
essential  aspects  of  the  primary  zones  of  turbojet 
combustors . 

2 .  RECOMMENDATIONS 

The  present  study  offers  the  following  recommendations  for 
further  activity.  These  complement  the  recommendations  offered 
in  Volume  I  of  this  report  and  in  References  2  and  3. 

•  The  Reynolds-averaged  computational  procedure  employed  in 
the  present  study  must  be  evaluated  in  a  number  of 
aspects.  These  include  grid  sensitivity  of  the  computed 
solutions,  effect  of  geometrical  and  fluid  mechanical 
parametric  variation  on  the  numerical  convergence 
(especially  in  reacting  flows),  and  the  influence  of 
variation  of  the  temperature  lower  bound  on  the 
predictions  of  reacting  flows. 

•  Some  of  the  interesting  implications  of  the  present 
predictions  and  their  flowfleld  ramifications  must  be 
investigated  through  an  integrated  experimental  and 
predictive  research.  Of  particular  significance  are  the 
CFD  examination  and  the  experimental  verification  of  (a) 
the  influence  of  bluff-body  blockage  ratio  variation;  (b) 
the  effect  of  the  variation  in  the  inner  and  outer  radii 
of  the  swlrler;  and  (c)  the  effect  of  variation  in  the 
swirl  angle. 

•  The  Reynolds-averaged  formulation  for  the  swirling  flows 
must  be  examined  in  terms  of  full  three-dimensional 
equations.  In  other  words,  the  validity  of  assumption  of 
axlsymmetry  in  the  mean  used  in  the  present  predictions 
remains  to  be  verified. 


It  is  likely  that  even  the  mean-field  predictions  are 
expected  to  suffer  from  the  anisotropic  effects  of  the 
near  wake  which  are  not  taken  into  account  in  the  k-e 
model.  A  first  step  to  consider  with  the  full  three- 
dimensional  formulation  is  the  incorporation  of  the 
algebraic  stress  model  (ASM) .  With  the  validated 
experience  of  ASM,  future  efforts  must  explore  the 
consideration  of  the  full  Reynolds-stress  equations. 

There  is  reason  to  believe,  unfortunately,  that  there  is 
no  guarantee  of  thermodynamic  realizability  connected 
with  the  stress-equation  model.  Nonetheless,  this 
question  remains  to  be  addressed  by  a  CFD  Investigation. 

Computational  aspects  relating  to  the  reacting  flowfields 
in  the  existing  procedure  (Reference  5)  remain  to  be 
fully  explored.  Thus,  even  with  the  k-e  model,  further 
research  is  needed  in  the  proper  Implementation  of  (a) 
the  Magnuasen-Hjertager  combustion  model,  (b)  the 
hydrodynamics-chemistry  coupling,  and  (c)  more  realistic 
chemical  kinetics. 

Finally,  the  proposed  intrusive  and  nonintruslve 
diagnostic  measurements  in  the  new  POSF  configuration 
must  carefully  examine  the  flowfield  implications  of  the 
present  predictions  to  help  arrive  at  the  modus  operand! 
of  experimental  verification  and  code  validation. 

Although  this  recommendation  is  addressed  to  the 
experimentalist,  it  has  a  significant  future  impact  on 
the  choice  of,  and  refinements  in,  the  several  CFD 
avenues  outlined  in  the  earlier  recommendations. 
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